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Abstract 

The equivalence is found between high-energy QCD in the generahzed leading loga- 
rithmic approximation and the one-dimensional Heisenberg magnet. According to Regge 
theory, the high energy asymptotics of hadronic scattering amplitudes are related to sin- 
gularities of partial waves in the complex angular momentum plane. In QCD, the partial 
waves are determined by nontrivial two-dimensional dynamics of the transverse gluonic 
degrees of freedom. The "bare" gluons interact with each other to form a collective 
excitation, the Reggeon. The partial waves of the scattering amplitude satisfy the Bethe- 
Salpeter equation whose solutions describe the color singlet compound states of Reggeons 
- Pomeron, Odderon and higher Reggeon states. We show that the QCD Hamiltonian for 
reggeized gluons coincides in the multi-color limit with the Hamiltonian of XXX Heisen- 
berg magnet for spin s = and spin operators the generators of the conformal SL{2,<C) 
group. As a result, the Schrodinger equation for the compound states of Reggeons has a 
sufficient number of conservation laws to be completely integrable. A generalized Bethe 
ansatz is developed for the diagonalization of the QCD Hamiltonian and for the calcula- 
tion of hadron-hadron scattering. Using the Bethe Ansatz solution of high-energy QCD 
we investigate the properties of the Reggeon compound states which govern the Regge be- 
havior of the total hadron-hadron cross sections and the small— x behavior of the structure 
functions of deep inelastic scattering. 
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1. Introduction 



The Regge problem is one of the longstanding problems in high-energy physics which goes back 
to sixties In a general sense, it is the problem of behavior of hadronic scattering amplitudes 
in the limit of high-energies s and fixed transferred momentum t 

s:^ -t^ M"^ (1.1) 

with M being a hadronic mass scale. There are a wealth of experimental data which still con- 
firm its intricacies, e.g., remarkably linear Regge trajectories of hadronic families, approximate 
duality of hadronic scattering amplitudes |Q, the growth of total hadronic cross sections and 
the small— a; behavior of the structure functions of deep inelastic scattering 0. 

There were numerous attempts to solve the Regge problems during the pre-QCD era which 
have led to the development of the theory of the Regge poles |1|]. The theory was inspired 
by the remarkable observation by Regge that, in quantum mechanics, the high-energy asymp- 
totics of the scattering amplitudes in potential scattering is governed by the singularities of 
the partial waves in the complex angular momentum plane. It was understood from general 
principals of relativistic quantum field theory [0] that a similar phenomenon occurs in high- 
energy hadronic scattering in Regge kinematics ( |1 . 1|) . Namely, after the decomposition of the 
scattering amplitude A{s, t) over partial waves as 

the behavior of A(s, t) in the Regge limit ( |1.1| ) depends on the singularities of the partial waves, 
A{uo,t), which lie in the complex l<j— plane to the left from the integration contour. Possible 
singularities of A{(jj,t) were identified as Regge poles and Regge cuts. Among all families 
of Regge poles there is a special one corresponding to the so called Pomeranchuk poles or 
Pomerons, which have vacuum quantum numbers and which provide dominant contributions to 
the hadronic scattering amplitude. However, Regge theory itself does not allow us to calculate 
the positions of these singularities and it is still a challenge for QCD as a theory of strong 
interaction to give us this input from "first principals" . After almost twenty years of a great 
deal of activity, we have different phenomenological models to describe the Regge behavior of 
the scattering amplitudes ||, ^, ^, but we do not have a complete understanding from QCD. 

In this paper we concentrate on perturbative QCD approach to the Pomeron [§. Namely, we 
consider hadrons as consisting of partons whose distributions are described by the nonperturba- 
tive hadronic wave functions in the infinite momentum frame. Hadrons scatter by exchanging 
quarks and gluons in the t— channel and we treat their interaction using the S*— matrix of per- 
turbative QCD. The hadronic scattering amplitude is a convolution of hadronic wave functions 
and partonic scattering amplitude. Although it cannot be calculated entirely perturbatively, 
we can find its evolution of as a function of energy s in perturbative QCD. The corresponding 
equations are called the evolution equations. 

There are at least two reasons why perturbative approach might make sense. First, since we 
do not understand dynamics in the nonperturbative regime of QCD, the perturbative approach 
is a first approximation to start with. Second, there is a simple way that perturbation theory 
itself could tell us about nonperturbative effects. Namely, studying ambiguities of perturbative 
series associated with the contribution of the so-called infrared renormalons one can get some 



insight into the structure of nonperturbative corrections 0, |T0|, [TT], |l^ . 
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Analyzing the scattering amplitude in the Regge limit ( p. . 1|) , one might be especially inter- 
ested to consider the case t = 0, corresponding to the elastic forward scattering. Then, the 
optical theorem |^ allows us to relate the imaginary part of the scattering amplitude to the 
total hadronic cross section 

<X^^) = ;^"^^V'^b(^'^ = 0)' (1-3) 



with Ah^h-g the amplitude of elastic scattering of hadrons /i^ and h^. Moreover, replacing one 
of the incoming hadrons by virtual photon with = —Q^ we get an analogous relation^ for 
the structure function of deep inelastic scattering 7* + /i — X 

F{x,Q^) = -lmA^*h{s,t = 0), s = Q^x , (1.4) 
s 



where —Q^ is a photon virtuality and x is the Bjorken variable. The Regge limit ( |1 . Ij) corre- 
sponds to s or, equivalently, to small— x asymptotics of the structure function F{x, Q^). 
We stress that relations ( |l.tj| ) and ( |1.4| ) have been found under the additional conditions of 



unitarity and analyticity of the scattering amplitude, A(s, t), in its different channels and that 
one of the consequences of these conditions is the Froissart bound 

crhth2i^) - const, log^ s , F{x, Q"^) < const, log^ x , s ^ 00 , x ^ . (1.5) 



Thus, in order to apply ( |L3D and (|1.4|) in perturbative QCD, we have to calculate hadronic 
scattering amplitude and preserve the conditions of unitarity and analyticity. 

Expanded in powers of the coupling constant as, the 2^2 scattering amplitudes have the 



following general form ||T^, |TB| in Regge kinematics ( |1 . 1| ) 



00 



A{s,t) = 2 5] «ss [{aJogsy^fmAt) + as{aJogs)"'-'fm,m-iit) + ... + aTfm,o{t)\ + 0(s°) , 

m=0 

(1.6) 

where the fm,i{t) are some functions of the transferred momentum. In the case of deep inelastic 
scattering, the C(s°) term includes logQ^ corrections to the scattering amplitude A^*h- Since 
we are not able to resum corrections to all orders in as, we need an effective approximation 
to the perturbative series which, first, will correctly describe the high-energy asymptotics of 
A{s, t) in perturbative QCD and, second, will be consistent with the unitarity and analyticity 
of the 5"— matrix. Combined together, these two conditions put severe restrictions on possible 
approximate schemes M, IT5|, One can easily satisfy the first condition by taking the 



leading logarithmic approximation 

a^logs/M^ ~ 1 , as<l, (1-7) 

in which only the fm,mit) term is kept in ( |1.6| ) and all nonleading terms fm,m-i{t), ■■■ , fmflit) 
are neglected. However, the resummation of leading logarithmic terms leads to the result [|l7l 



^lla(s, t)=if^ «,s(«, log s)™/„,„(t) (X (1.^ 



m=0 



^The exact relation between electron-proton inclusive cross section and the proton structure function 
F2{x, Q"^) can be found in 
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with Abfkl(O) = log 2 > 0, which violates the unitarity bound (|1.5|) . The s— channel uni- 
tarity of the S"— matrix is explicitly broken in the leading logarithmic approximation [|17|. To 
restore unitarity one has to take into account nonleading logarithmic corrections to the scatter- 
ing amplitude (pTol) which will involve fm,m-iit), fm,m-2(t), ■■■ terms. Thus, to understand the 
Regge behavior in perturbative QCD we need a new approach beyond the leading logarithmic 
approximation, which will preserve the unitarity of the S"— matrix. It is this requirement which 
makes the Regge problem so complicated in perturbative QCD. 

One of the possible ways to cure the unitarity problem, the so-called generalized leading 
logarithmic approximation, has been proposed in [0, |TB[ and can be formulated as followsQ 
Although the leading logarithmic terms fm,m(t) violate unitarity of (p!^), the total sum of 
both leading and nonleading terms satisfies the unitarity conditions. Therefore, there should 
exist some minimal number of nonleading terms in ( |1.6| ) which compensate nonunitarity of 
the leading logarithmic result. The generalized leading logarithmic approximation allows us to 
identify the "minimal" subset of such nonleading terms /mm~25 ••• • Of course, this 

procedure gives us only some pieces of nonleading terms fm,m-i, fm,m-2, ■■■ and not their exact 
expressions. The question remains on how small are nonleading corrections, / — /™™, which 
are left after a minimal set has been found. We have an answer |jl6| confirmed by lowest order 
perturbative calculations but the general proof is still missing. 

In what follows, we use the generalized leading logarithmic approximation to understand 
the Regge behavior of the scattering amplitude. We show that in this approximation QCD 
turns out to be closely related to two-dimensional exactly solvable spin chain models [1^, 



Such kind of relation has been expected for a variety of reasons pO| , pi] , p^ , [23| , P). The 
first calculations to lowest orders of perturbative QCD showed that the contributions of 
individual Feynman diagrams to the partonic scattering amplitude (p..6|) in the Regge limit 



( |1.1| ) have a form of effective two-dimensional amplitudes in which as log s plays a role of 
an expansion parameter and the dependence of fm,m on the transferred momentum t comes 
from two-dimensional integrals over transverse components (a = 1, 2) of the momenta 
of exchanged gluons.0 Another argument came from the physical picture of hadron-hadron 
scattering , which suggests that the Regge asymptotics is determined by a two-dimensional 
dynamics of transverse partonic degrees of freedom. In this paper we confirm these observations 
and identify the two-dimensional structure behind the high-energy QCD in the Regge limit as 
corresponding to the celebrated one- dimensional XXX Heisenberg magnet for the special value 
of the spin s = corresponding to a unitary representation of the conformal SL{2,€) group 



T9|. This correspondence allows us to apply powerful quantum inverse scattering method 



24| , |25| , P6i |2^ for the calculation of the nonleading logarithmic corrections /™° 



m,m— 2' ■■■ 

and it opens the possibility of finding the Regge asymptotics of scattering amplitude in 
QCD by means of the generalized Bethe Ansatz. 

The paper is organized as follows. In Sect. 2 we describe the generalized leading logarithmic 
approximation and find expressions for large— s behavior of hadronic scattering amplitudes and 
small— X asymptotics of the structure functions of deep inelastic scattering in terms of wave 
functions and energies of the compound Reggeon states. In Sect. 3 we review the properties 



^We should note that the generahzed leading logarithmic approximation restores unitarity of the scattering 
amplitude in its direct s, t and u— channels but not in the subchannels corresponding to the intermediate states 
0- 

•^Here we refer to the Sudakov decomposition of the gluon momenta over longitudinal and transverse mo- 
menta, K — aPj^ + l3Pg + k with and P^ hadron momenta. 
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of the Reggeon interaction and consider simplifications which occur in the hmit of multi-color 
QCD. In Sect. 4 we show that the holomorphic and antiholomorphic Reggeon hamiltonians 
coincide with the hamiltonian of XXX Heisenberg magnet for noncompact spin s = 0. The 
Schrodinger equation for the compound Reggeon states turns out to be completely integrable 
and the quantum inverse scattering method is applied to identify hidden conservation laws. In 
Sect. 5 we describe the generalized Bethe Ansatz which allows us to find the spectrum of the 
Reggeon states using the solutions of the Baxter equation. In Sect. 6 we discuss solution of 
the Baxter equation and general properties of the Reggeon states. Sect. 7 contains concluding 
remarks. 



2. Generalized leading logarithmic approximation 

We begin the construction of the scattering amplitude A{s,t) in the generalized leading loga- 
rithmic approximation by summarizing the properties of the leading logarithmic approximation 
( |1.7| ) in which only the coefficient functions fm,m(t) survive in (|T]^). 

In the leading logarithmic approximation [|, |T^, the partonic scattering amplitude is dom- 
inated by the contribution of soft gluons propagating in the t— channel between hadrons and 
interacting with each other. This interaction leads to the remarkable property of gluon reggeiza- 



tion IP, [13. The "bare" gluons interact with each other to form a collective excitation, the 
Reggeon, which has the quantum numbers of the gluon. However, in contrast with the gluon, 
the properties of the Reggeon depend not only on its momentum but also on the energy of 
particles interacting with the Reggeon. Introducing the Reggeon as a new elementary excita- 
tion, it is natural to replace QCD in the leading logarithmic approximation by an effective field 
theory in which the calculation of the scattering amplitudes can be performed in terms of the 
propagation of Reggeons in the t— channel and their interaction with each other. Each diagram 
in the effective theory is equivalent to an infinite sum of Feynman diagrams involving bare 
gluons. Then, the asymptotic behavior of the scattering amplitude is given in the leading log- 
arithmic approximation by the contribution of diagrams with only two Reggeons propagating 
in the t— channel, the famous ladder diagrams p, p^ . 

However, two Reggeon diagrams do not satisfy the s— channel unitarity condition and in the 
generalized leading logarithmic approximation one has to add additional diagrams to restore 



unitarity 0, |T5|, jTB] . This minimal set of diagrams is obtained from the two Reggeon diagram 



by iterating the number of Reggeons in the t— channel. For example, the first nonleading 
correction corresponds to the diagram with three Reggeons propagating in the t— channel. 
Continuing this procedure, we come to the conclusion that the scattering amplitude is given in 
the generalized leading logarithmic approximation by the sum of the diagrams shown in fig. |l[ 
These diagrams have a form of generalized ladder diagrams M, ITB|, in which summation is 



performed over all possible numbers of rungs representing the Reggeon interaction, and over all 
possible number of Reggeons in the t— channel, n = 2, 3, ... . It will be shown below that the 
scattering of each two Reggeons is described by the same effective theory as we started with in 
the leading logarithmic approximation. 

One of peculiar features of the Reggeon scattering is that it does not change the number of 
Reggeons in the t— channel. For n = 2 the diagram of fig.0 represents the leading logarithmic 
result for the scattering amplitude and it contributes to the coefficient functions fm,m{t) to all 
orders of perturbation theory. The contribution of the diagram with n = 3 is suppressed by 
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a power of as with respect to that for n = 2 and it determines the first nonleading coefficient 
/m,m-i ( |1-6D . In general, the contribution of the diagrams of fig.|I| to the partial wave A{uj, t) 
can be represented in the following form 

oo 

A{uj, t) = A^icu, t) + aA^ito, t) + ... = J2 ar'^n(c^, t) . (2.1) 

n=2 

Here, the n— th term corresponds to the diagram with n Reggeons in the t— channel, and the 
functions An{uj,t) have the following expansion 

A4u,t) = Y.(^) n fn-2Mi) ■ (2.2) 

Each term in this sum has a singularity at = which after integration in (|1.2|) is transformed 
into as(«slogs)' correction. When (|2.2|) and ( |2.1| ) are substituted into ( |1.2|) , they reproduce 



the perturbative expansion (|1.6|) . It does not mean however that the total sum, An{uj,t), has 
singularities only at = 0. If the partial wave An{uj,t) sums up to a singularity at u; = ujq, its 
contribution to the scattering amplitude (|1.2| ) would have been a Regge asymptotic behavior 
~ s'^°. Indeed, as we shall see below, this turns out to be the case for A2{uj,t), which has a 
square- root singularity at u; = Abfkl = log 2, the so-called BFKL Pomeron [|17|. We have 
to generalize this result for the higher partial waves An{uj,t) with arbitrary n and find their 
singularities in the complex cj— plane. 

2.1. The Bethe— Salpeter equation for the partial waves 

The partial wave An{uj,t) can be represented in the following form p, in terms of partonic 
distributions inside hadrons, $^ and $3, and partonic scattering amplitude, T„, 

A4u,t) = j[(fk] j[(fk'] $A({A:};g) Tn{{k},{k'}-u) ^^{{k'}-q), (2.3) 

where [(Pk] = Y[j=i (PkjS{J2j kj — q) and the definition of [d?k'] is analogous. Here, the hadrons 
hj^ and h-^ emit n Reggeons in the t— channel with momenta {k} = {ki, . . . , kn) and {k'} = 
{k[, . . . , k'^), respectively, and the total transferred momentum is g = Y^]=i kj = Sj=i with 
q^ = —t. The functions $a and $b describe the internal structure of the hadrons, and they 
can be expressed in terms of hadronic wave functions in the infinite momentum frame. In 
the high-energy limit (|1 . 1|) , both functions depend on the transverse momenta of Reggeons 



P, |T7| and not on the energy s, or equivalently on lu. The function T„ describes n to n 
scattering of Reggeons in the t— channel, and is the main object of our consideration. In 
(|2.3| ), the integration is performed over transverse momenta of n Reggeons, while integration 
over longitudinal components is performed inside T„(a;). The partonic distributions and the 
Reggeon scattering amplitude depend on the color indices of Reggeons and summation over 
these indices is implied in (^.3]). 

Considering the partonic distribution functions as states in a two-dimensional transverse 
phase space for the n Reggeons, one can rewrite the scattering amplitude ( |2.3| ) as the following 
matrix element 

Aniu^t) = ($A(g)|i;(a;)|$B(g)), t = -q\ (2.4) 
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where the transition operator Tn{uj) describes the elastic scattering of n Reggeons. To find the 
transition operator T„ we notice that diagrams of fig.|] have a ladder structure, which suggests 
a Bethe-Salpeter-like equation for Tn{uj) shown in fig.^ The corresponding equation has the 



following form |14, 28 



^T„M=tW+7^„T„H, (2.5) 

where T^^-* corresponds to the free propagation of n Reggeons in the t— channel, while the oper- 
ator Hn describes the pair-wise interactions of n Reggeons. The Reggeon scattering amplitude, 
Tn{{k}, {k'};uj), depends on the transverse momenta and the color indices of n incoming and 
n outgoing Reggeons. The hamiltonian Tin acts in ( p.5| ) only on the incoming Reggeons, and 
is given by 

^n = -^ E H^^nth (2-6) 

'^'^ n>i>j>l 

where the sum goes over all possible pairs of Reggeons. Each term in this sum has a 
color factor, which is given by the direct product of the gauge group generators in the adjoint 
representation of the SU{N) group, acting in the color space of i—th and j— th Reggeons, 



abc 



with fabc the structure constants of the SU{N). The Reggeon interaction ( p.6|) is described by 
the two-particle hamiltonian Hij, which depends only on the transverse momenta of Reggeons. 
If {/ci, /c„} and {k[, k'^} are the transverse momenta of incoming and outgoing Reggeons, 
respectively, then the operator Hij acts on Reggeon momenta as follows 



{ki, ...,kn\Hij\k[, ...,k'^) = H{ki,kj\k'i,k'j) Y[ 6'^{ki-k'i), (2.7) 



where H{ki,kj\k[^ kj) is given by |^ 

H-^H(tk-\k' k') - ^i mp-kT + ip-hn--ih-kdV .28) 



- 5''{ki - k[) J d^k' 



k''^ + {ki-k'f {p - k'f + {ki - k'f 



with p = kj + ki = k'j + k'^. The operator T^°) is equal to the product of n Reggeon propagators, 
and in momentum representation it can be written as 

{k„ k^\Ti'^\k[, k'^) = n S\k, - k'j)^ . (2.9) 

Iteration of the equation (|2.5| ) reproduces the ladder in fig. |l], and gives the perturbative expan- 
sion of Tn{uj) as series in Os/uj which leads to ( |2.2| ). But this is not what we are interested in 
because we want to resum the whole series in ( p.2| ). Using the Bethe-Salpeter equation (|2.5|) , 
we find the general solution for the transition operator 

Tniu;) = ^^Ti'K (2.10) 
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We recall that the hamiltonian 7i„ corresponds to the pair-wise interaction (|2.6|) of n Reggeons 
and it describes the evolution of the n— Reggeon state in the t— channel. 

Suppose for a moment that we know the spectrum of the Reggeon hamiltonian 

n })=E n,{a} \ Xn, W), (2.11) 

with {a} some set of quantum numbers, which parameterize possible solutions of the equation. 
Then, the eigenstates of Hn can be identified as compound states of n Reggeons, \xn), and the 
corresponding eigenvalues. En, determine the energies of these states. A simplest example of 
such states to be rederived below is the BFKL Pomeron which is built from n = 2 Reggeons. 



Once we solve (|2.11|) , we can expand the transition amplitude ( p.lO|) over eigenstates of Reggeon 
hamiltonian a^ 

^«M = E 1 IXn,w)(Xn,w|ri°), (2.12) 

where the sum over a means the summation over discrete and integration over continuous 
a. Combined with ( ^2.41 ), this expression implies that the singularities of the partial waves 
An{uj,t) in the cj— plane are determined by the eigenvalues En.{a} of the Reggeon hamiltonian 
Hn- Moreover, substituting (|2.12|) into (|2.4| ) and ( |1.2D , we formally take a residue at u; = Hn 
and get the following expression for the scattering amplitude 



oo 



Ais,t) = J2a:~'Anis,t), 

n=2 

An{s,t) = ^s($a|s^"TW|$b) = tsY.A"Ht)(3it\t) s^-w , (2.13) 

{a} 

where the residue functions and measure the overlapping between partonic distribu- 
tion functions and the wave functions of the compound states of Reggeons 

f3i^\t) = A(g)|Xn,w) , Pi'-Ht) = (Xn,w|T(°)|$B(g)) , it = -q') . (2.14) 

Although the definitions of the residue functions look different, they are closely related to each 
other as we will show below in (|5.29|) . We recall that the scalar product of states, used in (|2.4|) , 
( |2.13| ) and ( |2.14| ), implies integration over transverse momenta ki, ... , kn of n Reggeons as well 



as summation over their color indices. 

Thus, the scattering amphtude ( |2.13| ) has Regge behavior, An{s,t) ~ with the inter- 

cept related to the energy of compound states. Moreover, the scattering amplitude An{s,t) 
gets its dominant contribution from the compound states with maximum energy E^'^^ = 
ma.X{a} En^{a}- Ouce we have expressions for the wave functions of the Reggeon states we can 
use either phenomenological or model hadronic distributions to evaluate f3l^\t) and (3j^\t). 

The scattering amplitude ( p.l3| ) is given by a sum of An{s,t) over all possible numbers 



of Reggeons in the t— channel. With each individual An{s,t) we associate the family of the 
compound Reggeon states that appear as solutions of the eigenstate equation (|2.11|) . Thus, the 
original problem of calculating the scattering amplitude A{s, t) in the generalized leading log- 
arithmic approximation is reduced to the problem ( p.ll|) of the diagonalization of the Reggeon 
hamiltonian Tin for an arbitrary number of Reggeons. 
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To give a meaning to this expansion one has to specify the scalar product on the Hilbert space of the 



eigenstates and find the hcrmiticity properties of the hamiltonian. This will be done in sect. 5.3 
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2.2. Evolution equation for the structure function 

Let us apply ( p.l3| ) in the case of forward photon-proton scattering, A^*p{s,t = 0). Using (|T] 
and taking the imaginary part of the scattering amphtude ( |2.13| ), we obtain an expansion for 
the structure function F{x, Q"^) as a sum of nonleading corrections similar to that in (|2.13|) . 
The expressions for Fn follow from expressions ( p.l3| ) for An{s = Q'^/x,t = 0). Finally, we find 
the small— a: asymptotic behavior of the structure function of deep inelastic scattering in the 
generalized leading logarithmic approximation as0 



Fix, Q') = cyr'F^x, Q') + Oix' 



12\ _ „.n-2 T7 ( ^ , /-o^^O^ 

n=2 



F„(x,Q2) = E^"''"'*"*('^'r(Q^)(0)IXn,w) , (2.15) 

{a} 

where the Q^— dependence comes from $^.(q2)(0) and from 0{xP^ terms. To make the corre- 
spondence with the parton model, we represent Fn{x,Q'^) as convolution of a partonic cross 
section ctq and a partonic distribution function /, 

F„(x,Q2) = Jl^k] ao"-'^"(Q2;fci,...,fcO r--""(a;;fci,...,fc„) + 0(x°), (2.16) 

where [(fk] = (fki...(fknS^'^\ki + ... + kn), with a delta-function included to ensure the condition 
t = 0. The partonic cross section, ctq, describes the coupling of n Reggeons to the virtual photon 
through the quark loop. For large invariant mass, —Q^, it can be computed perturbatively as 
a scalar product of the photon state $^.((32)(t = 0) and the state of n Reggeons with transverse 
momenta ki, ...,kn and color indices 

a;ji-""(Q2; ki, fc„) = ($^*(Q2)(0)|/ci, ai; a„) . 

The partonic distribution, /, takes into account all effects of Reggeon interactions and the cou- 
pling of Reggeons to the nonperturbative hadronic state. In the generalized leading logarithmic 
approximation it is given by 

r'-''"{x;ku...,kn) = (A;i,ai;...;A;„,a„|a;-^"Tf)|$p(0)), (2.17) 

where fc„ = — (fci + ... + A;„_i) and the Bjorken variable x is small in the Regge limit. In terms 
of the eigenstates and eigenvalues, ( |2.11D , the distribution function can be written as 

r-'^"(x;fci,...,fc„) = 5:x-^-w x:Yo"}"(^i,-,U /?i"^(0), (2.18) 

{a} 

where we have used the following notation for the wave function of the compound state of n 
Reggeons, 

Xn^{a}"(^l5 •••5 ^n) = (^1; 0,1] kn, (ln\Xn,{a}) ; 

and /?^°^(0) is nonperturbative residue factor. Although we cannot compute the distribution 
function / in perturbative QCD, we may use expression ( |2.17| ) to find the evolution equation 
for / as a function of x 

^Y/n\^] m = -j [d'k'] n\^^'^'\{ky, {k']) /W(x; {k']) (2.19) 



^Here w e assu med that the energies of the compound states, En^{a}^ and the product of the residue functions 
are real in ( 2.13 ). The proof will be given in Sect. 5.3. 
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{ki, kn) and where the Reggeon hamiltonian was defined 



2 the evolution equation ( p. 191) coincides with the BFKL 



where {a} = (ai,...,a„) and {k} 
in (U), Q and For n 

equation |T7| . 

Comparing expressions ( p.l3|) and (p.l5|) , we conclude that the Regge asymptotics is defined 
in both cases by the properties of the Reggeon hamiltonian. The hadronic scattering amplitude 
and the structure function of deep inelastic scattering are given by an infinite sum of terms 
corresponding to the diagram of fig. |1| with fixed number n of Reggeons in the t— channel. The 
first term, n = 2, describes the leading logarithmic asymptotics and, taken alone, it violates the 
s— channel unitarity of the S"— matrix. Each next term in the sum over n in ( |2.13| ) and ( |2.15| ) 
defines a nonleading contribution, which is suppressed by a power of the coupling constant with 
respect to that of the (n — 1)— st term. Examining the high-energy behavior of the n—th term, 
An{s,t) and we recognize that the contribution of compound Reggeon states with 

positive energy > to (|2.13|) and (|2.15|) grows as a power of energy, and thus violates the 
Froissart bound ( |1.5| ). We recall however that these nonleading terms have been defined from 
the very beginning in such a way that to restore unitarity of the total scattering amplitude. 
This means, that, although each term in the sum ( ^.13 ) and ( ^.15| ) violates the unitarity bound, 
unitarity is restored in their sum.0 



3. Properties of the Reggeon hamiltonian 



We turn now to the solution of the eigenstate problem ( p.llD and (|2.6| ) in order to find the 



spectrum of the compound Reggeon states. The Reggeon hamiltonian (|2.6| ), ( p.7| ) and 
has been found from the analysis of Feynman diagrams contributing to the hadronic scattering 
amplitude in the generalized leading logarithmic approximation, and it inherits the properties 
of high-energy QCD in the Regge limit. 



The solution of ( |2.11| ) is known only for n = 2 and the corresponding compound state 



of two Reggeons is called the BFKL Pomeron |17|. The Reggeon hamiltonian was defined in 



(|2.6| ), (|2.7| ) and (|2.8| ) in the two-dimensional space of transverse momenta of Reggeons. It is 
more convenient however to analyse equation ( p.ll| ) in two-dimensional coordinate space, the so 
called impact parameter space, rather than in momentum space. To this end, we perform a two- 
dimensional Fourier transformation, and replace in (|2.3| ), ( |2.7| ) and ( p.8[ ) the two-dimensional 
transverse momenta ki,...,kn and k[,...,k'^ by two-dimensional impact vectors and 
b[, which describe the transverse coordinates of Reggeons. Then, for the impact vectors 
bj = {xj,yj) we define holomorphic and antiholomorphic complex coordinates ( 



z 



and analogous coordinates (zpz'j) for the impact vectors b'j. Now one can use ( p.7|) and 
to find the two-particle Reggeon kernel Hij in the impact parameter space. It turns out that, 
expressed in terms of holomorphic and antiholomorphic coordinates, the kernel Hij becomes 
holomorphically separable p9|, i.e.. 



Hij = H{zi, Zk) + H{zi, Zk) . (3.1) 



^To get some insight into possible mechanisms of unitarity restoration one may consider the following exam- 
ple, s'*' — + + ... = — A/aglogs, in which each individual term of the series violates the unitarity 
bound but their sum does not. 
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Here, two operators on the r.h.s. act separately on holomorphic and antiholomorphic coor- 
dinates of Reggeons. After Fourier transformation of ( |2.7D and (|2.8|) they are given by the 
following equivalent representations 

H{z„Zk) = -PrHog{zi-Zk)P^-Pk^\og{zi-Zk)Pk-\og{P,Pk)-2-^^ 

= -2\og{zi - Zk) - {zi - Zk) \og{PiPk){zi - Zk)'^ - 27jj , (3.2) 

where Pj = i-^ and 7^ is the Euler constant. The same operator can be represented as [|29[ 

°° 2/ + 1 2 d"^ 

Substituting this expression into ( p.lD and (p.6|) , we find that the Reggeon hamiltonian Hn is 
invariant under the conformal transformations ^ 

azi + b _ azi + b 

Zi^ T^' - , (3.4) 

czi + a czi + a 

with ad — be = ad — be = 1. Indeed, the generators of these transformations, 

n n n 

S' = Y.^kdk, S~ = -Y.dk, S+ = Y4dk (3.5) 

k=l k=l k=l 

and the analogous antiholomorphic generators S^, S~ and form the 5*1/(2, C) algebra and 
commute with L^^, and, as a consequence, with 7i„. 

Let us consider the properties of the eigenstate Xn,{a} of the Reggeon hamiltonian (p. 11 



m 



the impact parameter space. These states are parameterized by quantum numbers {a}, which 
should appear as eigenvalues of some "hidden" integrals of motion, and by a two-dimensional 
real vector bo, which represents the center of mass of the compound Reggeon state. In this no- 
tation, the wave function Xn,{a} = Xn,{a}{{bi}; &o) = Xn,{a}{{zi, Zi}; zq, zq) satisfies the relations 
(|2.11|), (|2.(j|) and (|3.1|), or equivalently 



En,{a}Xn,{a}{{Zi, Zi}] Zq, Zq) = ^ [H{Zj, Zk) + H{Zj, Zk)] tjtl Xn,{a}{{Zi, Zi}] Zq, Zq) , 

j,k=l,j>k 

(3.6) 

which can be interpreted as a two-dimensional Schrodinger equation for a system of n pair-wise 
interacting particles with coordinates {zi, Zi] and internal color degrees of freedom. One may 
try to rewrite the total hamiltonian in (|3.6|) as a sum of holomorphic and antiholomorphic 
parts using the fact that H{zj,Zk) and H{zj,Zk) commute. However, the resulting two terms 
do not commute with each other due to nontrivial color factor in ( |3.6|) , and as a consequence 
holomorphic and antiholomorphic sectors become coupled to each other via color degrees of 
freedom. It is this interaction which makes difficult the solution of the Schrodinger equation 

dp). 

Each Reggeon carries a color charge and the total charge of n Reggeons is equal to their 
sum J2]=it'j- Since the compound Reggeon states propagate in the t— channel between two 
hadrons, they carry zero color charge, unchanged by the Reggeon interaction, 

n n 
[^n,E^I]=0' E^"IXn,w) = (3.7) 

3=1 i=i 
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An essential simplification occurs in (|3.6|) in multi-color limit, ^ oo and asN = fixed, 
29| , |30| . In this limit, only planar diagrams of fig.|l| survive, which have the form of a cylinder 



attached by both edges to the hadronic states |3y] . Reggeons propagate along the sides of the 
cylinder and it makes them possible to interact only with two nearest Reggeons. Using the 
double line representation for Reggeon color charge and applying the standard rules of large 
counting, one finds that the color structure in ( p.6|) can be simplified asQ 

t-t^^^-N, forn = 2; t^t^ ^ -^Sij+i , forn > 3 , 

where i,j = 1, ...,n and the Reggeons with i = 1 and i = n + 1 are considering as coinciding. 
Then, in the multi-color limit we find that, first, the color factors become trivial in ( p.6| ) and 
as a consequence holomorphic and antiholomorphic sectors become decoupled and, second, 
inside each sector in ( |3.6| ) the pair-wise Reggeon interaction is replaced by a nearest-neighbour 
interaction with periodic boundary conditions. Thus, in the large— A^ limit, the two-dimensional 
Schrodinger equation (|3.6|) is replaced by a system of two one-dimensional Schrodinger equations 



29 



Hn^n,{a}{{Zi}; Zq) = £n,{a}^n,{a}{{Zi}; Zq) , Hn(Pn,{a}{{Zi}; Zq) = en,{a}Vn,{a}i{Zi}; Zq) . 

(3.8) 

The hamiltonians if„ and Hn are defined as 

n n 
Hn = ^H{Zk,Zk+l) , Hn = ^H{Zk,Zk+l) , (3.9) 

k=l k=l 



with two-particle hamiltonians given by ( |3.2D or ( PTB| ) and Zn+i = Z\. Once we know the solution 
of (|3.8|) , the eigenstates (|2.11|) of the Reggeon hamiltonian in the multi-color limit can be found 



as 



Xn,{a}{{Zi, Zi}; Zo, Zo) = (pn,{a}i{Zi}; Zq) y?„,|„} ({^i } ; ^o) , (3.10) 

and the corresponding eigenvalues are given by 

^n,{a} = -^;^[^n,{a} + £n,{a} J ■ (3-11) 

We identify (|3.8D as a system of Schrodinger equations for two one-dimensional lattice mod- 



els with nearest neighbour interaction (3^), and with the number of lattice sites equal to the 
number of Reggeons. The hamiltonians of the models are defined by holomorphic and anti- 
holomorphic parts of the Reggeon hamiltonian ( p.l|) , and the quantum space for each site of 
the lattice, h^, {k = l,...,n), is parameterized by holomorphic and antiholomorphic coordi- 
nates, respectively. Finally, we notice that and Hn are separately invariant under conformal 
transformations ( p.4|) . 

Thus, using the properties of the Reggeon interaction in the multi-color limit, we reduced 
the original problem ( p.ll| ) and (|3.6|) of the diagonalization of the Reggeon hamiltonian to 
the solution of the one-dimensional lattice models ( |3.8| ). However, this alone does not at all 
guarantee that these lattice models can be exactly solved. To explore this possibility we recall 



^There are two special case, n = 2 and n = 3, when the color factor can be calculated for finite N. For n — 2 



we use (3^) in the form ii + <2 = to get i"^2 = —titi = —Ca] for n — 3 a similar condition ti + 12 + ~ 
leads to t^q = l/2(<^ + q)'^ - l/2{tiq + qq) = -Ca/2 with Ca^ N the Casimir operator. 
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the famous example of the exactly solvable one-dimensional lattice model, the XXX Heisenberg 
chain of interacting spins s = 1/2 with the hamiltonian 



1 " 

iJ.=i/2 = -- E («+i - l) ' ^3.12) 
^ k=\ 

where o"^ are Pauli matrices acting in the fc— th site of the lattice. This model can be generalized 



JI| by means of the quantum inverse scattering method to give a family of XXX Heisenberg 
models for an arbitrary complex values of spin s. The unique feature of these models is that 
they are completely integrable, that is, that they contain additional integrals of motion whose 
number is equal to the number of degrees of freedom. The reason why we might be interested 
in considering complex spin s is that the one- dimensional lattice models (|3.8|) , (|3.9| ) and ( |3.3|) , 
which describe the Regge behavior of multi-color QCD, coincide with the XXX Heisenberg 
magnet for spin s = 0. Let us now discuss how we may make this identification. 



4. Multi-color QCD as XXX Heisenberg magnet 



To identify the lattice models ( |3.8| ) we forget for a moment about their QCD origin and construct 
the hamiltonian of the XXX Heisenberg magnet for spin s [^. Let us consider the one- 
dimensional lattice with n sites and assign the spin operators S"^ with a = 1,2,3 to all sites 



k 



l,...,n, 



Si 



4dk 



2szk 



-dk 



Zkdk 



(4.1) 



where S*^ = S*^ ± iS"^. The spins Sk act as differential operators on the local quantum space 
hk corresponding to the k—th site. The total spin of the lattice S = Y^^^i acts on the full 
quantum space H„ = (8>^=ihfc, and for s = the operators 5''* coincide with the generators of 
conformal transformations, eq. (|3.5| ). The definition of the integrable XXX spin chain is based 
on the existence of a fundamental operator Rfcm(A), which acts on the space hk^hm, and which 
satisfies the Yang-Baxter equation [0, |2^, |26|, p7| 

Rfcm(A - fi)Rkl{\ - p)Rm/(/" - P) = RmliP - p)RfcK-^ " p)Rfcm(A - /i) , (4.2) 

with A, fi and p arbitrary complex spectral parameters and k,m,l = 2+. The solution of this 
equation for an arbitrary complex s is given by [0, [3^ 



Rfcm(A) 



T{i\-2s)T{iX + 2s + l) 



r{iX - JkmWiiX + Jkm + 1) 

where the operator J^rn acts on the space ® and satisfies the equation 

"^fcm(l + Jkm) = {Sk + SmY = '^SkS-m + 2s(s + 1) . 

The hamiltonian of the XXX Heisenberg magnet for spin s is defined as 



(4.3) 



E Hk,k+1 , 
k=l 



Hk,k+1 = —'i'-TZ^^SRk,k+liX) 

dX 



(4.4) 



A=0 
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In the special case, s = 1/2, one can recover the hamiltonian (|3.12|) from this expression while 
the limit s ^ oo is related to the nonlinear Schrodinger equation. Let us consider the case 
s = 0. Using the explicit expressions ( |4.1| ) for the spin operators for s = we obtain the 
two-particle holomorphic hamiltonian up to an inessential (infinite) constant^ as 



k,k+l 



with ip{x) 



-■ip{-Jk,k+i) - V^(l + Jk,k+i) + 2^/'(l) , Jfc,fc+i(l + Jk,k+i) = -{zk - Zk+ifdkdk+i 

(4.5) 

we find that both expressions for the two- 



£r(a;) 



Comparing ( [4.5|) and (|3] 



particle hamiltonians coincide after we identify Jk,k+i{)- + Jk,k+i) = ^Ik+i perform sum- 
mation over / in (|3.3|) . This means |T9[ that the holomorphic and antiholomorphic Reggeon 
hamiltonians ( |3.9| ) are identical to the hamiltonian ( |4.5| ) of the XXX Heisenberg magnet for spin 
s = 0. Among other things, this remarkable property implies that the system of Schrodinger 
equations (|3.8| ), which describes the Regge asymptotics of multi-color QCD in the general- 
ized leading logarithmic approximation, is completely integrable, and that one can apply the 
quantum inverse scattering method to identify the "hidden" integrals of motion. 

To this end, we follow the standard procedure |2i, 25, ^ 27| and define the auxiliary and 
fundamental Lax operators, Lfc,a(A) and Lfcj(A), respectively, for all sites of the lattice. 



Lk,aW = A4 ® 4 + iSk ® <?a 



A + iSl iS^ 



A 



iSl 



= Rfc,/(A) . (4.6) 



These operators act in the space \ik®V with the auxiliary space V = 'C for Lk,a and V = )if 
for Lfcj (the auxiliary space hf has the dimension of the local quantum space h^,). The Lax 
operators satisfy the Yang-Baxter equations similar to that in (|4.2|) . Now, we multiply Lax 
operators as matrices in the auxiliary space and define the auxiliary monodromy matrix 



(A) = L„,,(A)L 

n~l,a 

(A)...Li,,(A) 



( A{X) 




and analogously the fundamental monodromy matrix 

TXA) = L„j(A)L„_ij(A)...Lij(A). 



(4.7) 



(4i 



Finally, we take a trace of the monodromy matrices over auxiliary space and obtain the auxiliary 
and fundamental transfer matrices 



A(A) = tr,T,(A) 



r(A) = tr;T;(A) , 



(4.9) 



which are operators in the full quantum space H„ of the lattice. Due to the properties of the 
Lax operators (^4.6| ), the operators r and A thus defined commute with each other for different 
values of the spectral parameters 



[r(A),r(/i)] = [r(A),A(/.)] = [A(A), A(/i)] = , 

From this relation we find that the system of Schrodinger equations (|3] 
mutually commuting local integrals of motions 



(4.10) 

has two families of 



Ih = —i 



logr(A) 



Qn—k 



A=0 



1 d'' 

~k\d)^ 



A(A) 



(A; =1,2, 



(4.1i; 



A=0 



®As we will show in Sect. 5.3 infinite constants c orresp onding to holomorphic and antiholomorphic hamilto- 
nians, Hk^k+i and Hk,k+i, are canceled in the sum 



14 



We conclude from (14TT|) and ( ^lOD that || 



[Hn,Tk] = [Hn,qk] = [qk,qj] = [Ik,Ij] = [qk,Ij] = 0. 



(4.12) 



Here, the operators Ik and q^ describe the interaction between {k + 1) nearest neighbors, and 
in the special case k = 1 the operator Xi coincides with the holomorphic Reggeon hamiltonian 
(I3T 



II — Hn — ^ Hk,k+1 ■ 



(4.13) 



k=l 



The expressions for Xk with k > 2 are more complicated, but one can easily find expressions 
for the operators qk by using the explicit form |l^, ^ of the auxiliary Lax operator ( ^^) for 
spin s = 0, 

Lk,a = XI + ®{zk,-l)dk. (4.14) 



Then, using the definition ( ^.101 ) we find the auxiliary transfer matrix as 

A(A) = 2A" + gaA"-' + q^X^-^ + . . . + g„ 
with the operators qk given by 



n>ii>i2>--->jfc>l 



with k = 2, ...,n and Zjk = Zj — z^. In particular, for q2 we have expression 



q2 



J2 z%djdk = -S^S"^ - ^{S+S- + S~S+) = -h{h - 1) 



(4.15) 



(4.16) 



(4.17) 



n>j>fc>l 



which we identify as the Casimir operator of the conformal group ( |3.5| ) with h being the 
conformal weight. Then, the relation [g2, H-n\ = follows from the invariance of the holomorphic 
hamiltonian under conformal transformations (|3.4| ). 

Thus, the holomorphic Schrodinger equations ( p.8| ) has a sufficient number of integrals of 
motion ([4.11|) and ( [4.12| ) to be exactly solvable. It is clear that the same consideration can 
be performed for the antiholomorphic sector in ( |3.8| ). The wave functions of the compound 
Reggeon state, Xn{a}, depend on the quantum numbers {a}, which we now identify as the 
eigenvalues of the integrals of motion ( [4.11| ) and their antiholomorphic partners. Moreover, the 
original complicated problem of solving the Schrodinger equations (p.8|) can be reduced to the 
diagonalization of the conservation laws ( [4.11[ ). We remember that in the case of the XXX 
Heisenberg magnet for spin s = 1/2 this has been done by means of Algebraic Bethe Ansatz 



2^ , |25| , |2^ , ^ . However, this method requires the existence of a highest weight for each site 
of the lattice. This makes it inapplicable for spin s = 0, corresponding to infinite-dimensional 
representations of the S'L(2,C) group. To solve the XXX Heisenberg magnet for spin s = 0, we 
have to generalize the Bethe Ansatz to the case of the noncompact conformal group. 



5. Generalized Bethe Ansatz for Reggeon Hamiltonian 

The Generalized Bethe Ansatz for the XXX Heisenberg chain for spin s = has been developed 
in Ref. |jT9|. It is based on the method of the Q— operator |3^, |3^ and on the relation 
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between the XXX magnets for spins s = and s = —1. To understand the latter relation we 
consider expression ( [4. 1]) for the spin operators for s = and s = — 1 and find the following 
identity, Sjf~ = Pk S'f.'*~°'*P^^ with Pk = idk- It leads to the similar relation between Lax 
operators 



4r'^(A) = p. Lt:'\\) p,' , Li:7"^)(A) = i.tf\\) p,-' , (s.i) 

where the additional factor in the last relation follows from the numerator of the fundamental 



P— operator ( |4.3| ). After its substitution into ( f4.7|) and (^4.8|) we find the relation between 
transfer matrices (|4.9|) in both models 



A(— 1)(A) = Pi...P, A(^=0)(A) (Pi...P„)-' = (Zi2^23-^nl)~'A(^=°nA) Z,2Z23...Znl (5.2) 

and 

/ ■ \ _ 1 \ " 

-'^'=~'\X) = Pi...P„ r(^=°)(A)(Pi...P„)-^ = (zi2^23-^nl)~V(^=°)(A) z,,z,,...z^, , 



iA + 1 



where we took into account that operator g„ = zi2Z23---ZniPi...Pn, defined in ([4.16|) , commutes 
with the transfer matrices (|4.9|) . Hence, the eigenstates of the transfer matrices in the XXX 
models for spin s = and s = — 1 are related as 

^t~'\{^^h^0) = iZuZ23-Znir'^t'\Ui};^0) = g^^PlPs-Pn ^^^^ ({^.}; ^o) • (5.3) 

For corresponding eigenvalues of the hamiltonians and the integrals of motion we have from 
dg), (^) and (^ 

et'^=ef'-^n, qf^ = ql'^'^ (5.4) 

and 

7-(s=0) _ 7-(s=-l) 7-(«=0) _ 7-(s=-l) _ o„/ \fc/oL\i 



where and stand in i f).3\) and (|5.4| ) for eigenvalues rather than for operators. These 
relations establish the equivalence between XXX Heisenberg magnets for spins s = and 
s = -1. 

We recall that the holomorphic Reggeon hamiltonian ( p.8| ) coincides with the hamiltonian 



(|4.4| ) and ([4.5| ) of the XXX magnet for spin s = and therefore their spectrum are identical 

^n^^\{Zi}; ^o) = Vn,{a}{{Zi}; Zq) , 4'"°^ = ^n,{a} 

It is clear that the same correspondence takes place for the antiholomorphic Reggeon hamilto- 
nian, Hn, in ( p.8|) . The spectrum of the Reggeon hamiltonian depends on the set of quantum 
numbers with the conformal weight h being one of them. As it follows from the definition ( ^.17| ), 
parameter {—h) has a meaning of the total spin of the XXX chain. In the case of the compact 



spin s = 1/2 XXX chain it can take either integer or half integer values [^, ^ |26|, For 
noncompact spin s = chain, possible values of h can be even complex and their explicit form 
will be found in Sect. 5.3. However, among all possible h there is a subset of integer positive 
h > n which plays a special role. 
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5.1. Algebraic Bethe Ansatz 

The reason why we prefer to deal with the spin s = — 1 in constructing the Bethe Ansatz is 
that there exists a pseudovacuum state \Q) in the total quantum space of the lattice which is 
the highest weight in each site 



^ 2 2 /y2 



sl\n) = -\n) 



with the spin operators given by ( |4.1| ) for s = —1. This property allows us to apply the Algebraic 

I, and construct the special subset of the eigenstates corresponding 



Bethe Ansatz m, |5l 

to the integer positive values of the conformal weight h > n in ( 4.17|) . The corresponding 



eigenstates, the so-called Bethe states, are defined as [|T9 



I = h — n = 0,1, 



(5.5) 



nO 



, Zn and 



Here, B(X) is the differential operator which acts on holomorphic coordinates zi 
which can be found through (|4.71 ) and using ( [4. 141) and (|5.1| ) as an element of the auxiliary 
monodromy matrix, T^(*"~^\ for the XXX magnet of spin s = — 1. As an example, we give 
expression for B{X) for n = 2 



Bn=2{\) = -A (Pi + P2) + P1P2 Z,2 . 

The Bethe states ( p. 5] ) are parameterized by the set of complex numbers {A} = (Ai, 
which are solutions of the Bethe equation for spin s = — 1 



A. 



Afe + i 



h—n 

n 



Xk - Xj + i 



X, 



A, 



(5.6) 

A/i-n) 

(5.7) 



The eigenvalues of the holomorphic Reggeon hamiltonian corresponding to the Bethe states 
( |5.5| ) are given by 

1 



h—n 



e 



The Bethe states 



a + 1 

diagonalize the transfer matrices 

h—n 



(5. 



\ _ \ _i_ n h-n \ 

A(A) = (A - iT n \ V + (A + ^r n \ — r 



A 



k=l 



and we can find the eigenvalues of the conserved charges {q} from this expression by expanding 
it in powers of A and comparing with ( ^.15| ). The eigenvalues of the operators {1} can be found 
from (|5.4|) using the relation 



h—n 



it-'^ = ^{-)\k - 1)! E [(A, + z)-^' - (A, - 



Given a solution of the Bethe equation ( |5.7| ), the relations ( p. 51 ) and ( |5.8| ) allow us to find the 



solution of the holomorphic Schrodinger equation in ( p.8|) , the holomorphic wave function (5.5) 
and the energy ( |5.8| ). We stress that these results correspond to the special case of positive 
integer conformal weight h > n while in general we are interesting to choose /i to be a complex. 
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5.2. Functional Bethe Ansatz 



The explicit form of the Bethe states ( ^.51) can be used in order to generahze them for an 
arbitrary complex h. The generalization is based on the existence Q of the Baxter operator 
(5n(A) which acts on the full quantum space of the XXX spin s = — 1 model and commutes with 
itself, [Qn(/^), Qn(A)] = 0, and with auxiliary transfer matrix, [Qn{fi), A{\)] = 0. Moreover, the 
eigenvalues of the operator satisfy the Baxter equation 



(2A" + gsA"-' + gsA"-^^ + ... + g„)Q„(A) = (A + ^)"Q.(A + ^ + (A - 0"g„(A 



(5.9) 



where qk and Qn{^) denote the eigenvalues of the corresponding operators in the XXX model 
for spin s = — 1. The relation ( p.4| ) allows us not to indicate explicitly the value of spin. Solving 
the Baxter equation we will get expression for (5„(A) which depends on the quantum numbers 
h and and which is the function of complex spectral parameter A. Then, for given solution 
Qn{^) of the Baxter equation the eigenvalues of the XXX magnet hamiltonian H^^~^^ and of 



the integral of motions 1^ 



{s=-l) 



are given by 



d\ ^Q„(-A + z) 



-1-1 



{s=-l) 



. , Q„(-A 

—t^—r log ■ 



A=0 



dA^ 



Qni-X + i) 



A=0 



Using the relation (|5.4|) between XXX models for spins s = and s 
of the holomorphic Reggeon hamiltonian 



-1 we find the spectrum 



= —1--r7 log — 

Qn(-A + z) 

and the eigenvalues of the integrals of motion 



Q„(A) = A"g„(A) 



(5.10) 



A=0 



^ .d\ Qn{-\-l) 

^k = -I— log 



d\^ 



Qn{-X 



A=0 



To define the corresponding eigenstates we introduce the Sklyanin operators for the XXX 
magnet of spin s = — 1. We use the definition (|4.6| ) and ([4.1|) of the auxiliary Lax operator for 
s = — 1 and find from ( [4.7| ) that the operator -B(A) is a polynomial of order n — 1 in A which 
we represent as 

B{X) = tS'iX - x,){\ - X2)...iX - x„_i) , (5.11) 



where S = — J2k=i 9k is the total spin of the model and xi, .. 
The relative order of factors is inessential in ( p.ll| ) since [S~ 



Xn-i are operator zeros of -B(A). 
a^fc] = [a;A:,a;j] = 0. The exphcit 



form of the operator zeros Xi, can be found for any fixed n from the definitions (|5.11 

and ( |4.7|) . In particular, using (|5.6| ) we find the explicit expression for the operator zero at 
n = 2 

^^^^-12. (5.12) 



Xi 



1 di + d2 

The operators -B(A) and Xk have different expressions for s = — 1 and s = models and the 
relation between them is similar to that between the auxiliary monodromy matrices, eq.( [5.2| ). 
Then, having expressions for the operators xi, in the XXX model for spin s = —1 we 



18 



use the solution of the Baxter equation ( pTOf ) to find the eigenstates of the holomorphic Reggeon 
hamiltonian in the following form 

Vn,{q}{{Zi}; Zo) = Zi2Z23...Znl {iS^Y'^'" (j^l )Qn (3^2) ■ ■ -Qn (^n-l) 2 2 2~ (5.13) 

^10^20 •••^raO 

with h a complex conformal weight. Notice that the operators Qn(A) and xj do not commute, 
and Qnixj) denotes the eigenvalue of the operator (5n(A) evaluated for operator value of the 
spectral parameter X = Xj. 

Expressions ( |5.10| ) and ( ^.131 ) give the solution of the holomorphic Schrodinger equation (|3.8| ) 



for an arbitrary complex conformal weight h. Different eigenstates of the Reggeon hamiltonian 
are parameterized by different solutions of the Baxter equation ( p.9| ). In particular, for integer 
positive conformal weight h > n the Baxter equation (|5.9|) has a polynomial solution 

h—n 

g„(A) = const, n(^-^'fc) (5-14) 

k=l 

and taking A = Xj in ( |5.9| ) we find that the parameters Ai, Xh-n satisfy the Bethe equation 
( |5.7| ). After substitution of ( |5.14| ) into ( |5.13| ) we use the definition ( |5.11| ) to identify ( [5.10| ) and 



(|5.13| ) with analogous expressions, (|5.8|) and ( |5.5|) , given by the algebraic Bethe ansatz. 



To apply the expressions ( |5.10|) and ( ^.13| ) one has to solve the Baxter equation (|5.9|) for an 



arbitrary values of the quantum numbers h and qk and under appropriate boundary conditions. 
In particular, the condition for the solution of the Baxter equation ( ^.91) to be a finite polynomial 
in A leads to the algebraic Bethe Ansatz described in the previous section. 

5.3. Scalar product 

The decomposition of the Reggeon hamiltonian over the complete set of the states has been 
introduced before in ( p.l2|) . However this decomposition is formal unless we specify the scalar 



product on the Hilbert space of the Reggeon hamiltonian. Although we have the general 
expression ( |5.13| ) for the compound Reggeon states, the choice of a scalar product together with 



the condition of finiteness of the norm may rule out some of these eigenstates and introduce 
selection rules for the quantum numbers h and qk- In what follow, we use the scalar product 



proposed in and based on the relation between XXX magnets for spin s = and s = —1. 



Let us define the following auxiliary hamiltonian 

n 

nr = Ht~'^ + Hf^ = E H^J+i , (5. 15) 

fc=i 

which is equal to the sum of holomorphic XXX hamiltonian for spin s = — 1 and antiholomorphic 
hamiltonian for spin s = 0. Using relations ( |5.4| ) we find the spectrum of this hamiltonian as 

paux _ (s=-l) _|_ =is=0) 

x'''''i{zk,Zk};zo,zo) = if^'=-'\{zk};zo)^^'='\{zk};zo). 

These expressions are related to the eigenvalues and eigenstates of the Reggeon hamiltonian, 
(^JOj ) and Q, as 

X''''''{{Zk,Zk}]Zo,Zo) = {zi2Z23--Zniy^Xi{Zk,Zk};Zo,Zo) . (5.16) 

= q~^PiP2--Pn xi{zk, Zk}; zo, zq) 
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Then, in each site of the model (|5.15| ) we have three holomorphic spin s = —1 operators 
and three antiholomorphic spin s = operators which form the unitary representation of the 
principal series, of the SL{2,€) group and we identify as the local quantum space 
hfc in the k—th site. This identification leads to the following properties. We notice that two- 
particle hamiltonian -ff^'^^,.! is unbounded operator due to singularities of the function on 
real negative axis. 



-^12 



-H-J12) - ^(1 + J12) - H-J12) - HI + J12) + c 



(5.17) 



where JuiJu + 1) = —did2zf2 and JuiJu + 1) = — ^i2<9i92 are two Casimir operators of 
the SL{2,€). The operators J12 and J12 act on the tensor product of two quantum spaces 
which can be decomposed into the direct sum of the principal series representations 

ger m and real u. For fixed m and u the Casimir operators have 



the following eigenvalues 



J 



12 



-h, 



12 



h* 



1 + m 



After their substitution into ( p.l7| ) we use the identity il){x) — '?/'(l — x) 
that singularities at = cancel in (|5.17 ) to give a finite real result 



rr aux 
-"12 



. ^ , /I + |m| ^ 
-AReipl — - + IV 



+ C 



-vrcotfvrx) and find 



(5.1^ 



The constant C entering into ( |5.17D is defined as| 



C = ij{-2i6) + + 2i6) + ij{2 - 2i6) + + 2i6) 



5^0 



2 + 4ij{l). 



(5.19) 



The total quantum space of the model is the tensor product of n copies of and the scalar 
product on this space is given by 



(5.20) 



where dzdz = YVk=i dzkdzk. With this choice of the scalar product the hamiltonian Ti^^ as well 
as the Reggeon hamiltonian. Tin, are bounded operators on the quantum space of the model. 



Using one-to-one correspondence, ( 5.16|) , between the spectra of the hamiltonians 7i^"^ and 



Tini we obtain two different choices for the scalar product of the Reggeon states x [HI 



j dzdz xiiz,z)PiPi...PnPniX2iz,z))* or J 



dzdz Xii^'_')iX2{z,_z)r 

Z\2Z\2---Zn,\Zn,\ 



For eigenstates of the Reggeon hamiltonian (|3.10|) both definitions are equivalent. 

Another important property of the hamiltonian ( |5.15| ) is that it is a selfadjoint operator. 



nja 



^To derive ( 5.17 ) and the expression for C one uses the definition (4^) and ( |4.3| ) of the two particle XXX 
hamiltonian for antiholomorphic spin s and for holomorphic spin (—1 — s*) with s = iS and 5 — > 
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since for s = — 1 and s = the spin operators (|4.1|) are related as S'"^|s=„i = 5'"|s=o- Together 
with ( gl6D this means that, first, the eigenvalues of the Reggeon hamiltonian corresponding to 
the eigenstates with the finite norm ( |5.20| ) are real 



Im E„ 



Air 



Im (sn + eJ 



(5.21) 



where En is the energy of the compound state of n Reggeons defined in ( p.ll|) and (|3.11 



Second, the eigenstates ( p.l(j| ) are orthogonal with respect to the scalar product ( |5.2(JD and, as 
a consequence, we have the following orthogonality condition 



Sq,q'6{zo - z'q)5{zq - z'q) = / dzdz 
and the completeness condition 



{Xn,{q'}{{Zk, Zk}', z'q, 2g)) Xn,{q}{{Zk, Zk}', Zq, Zq) 



Zl2Zl2Z23Z23---ZnlZnl 



(5.22) 



n „ 

5(2;^ - z'^.)5{zk - 4) = H / dzodzt 

k=l {q} 



Xn,{q}{{Zk, Zk}', Zq, Zq) (^Xn,{q}{{Zk, 4}? ^0; Zq)^ 



Zl2Zl2Z23Z23---ZnlZnl 



(5.23) 



where sum over {q} means summation over discrete and integration over continuum quantum 
numbers g^. 

We recall that relations ( |5.21| ) are valid only for the compound Reggeon states with the 
finite norm ( |5.20| ) and this condition restricts possible values of the quantum numbers {qk}- In 
particular, the energy of the compound state is given by ( p.ll|) 



with En = Sn{q) defined in ( |5.10|) and with analogous expression for = £„(?)• To find the 
energy one has to solve the Baxter equation ( |5.9| ) and it is natural to assume that for real 
{q} and A the solution of the Baxter equation, Q„(A), may be taken to be real. Then, it follows 
from ( p.lO[ ) and ( p. 91 ), that £«((?) is a function of holomorphic quantum numbers {qk} with real 
coefficients and, hence, (e„(g))* = £:„(g*). We notice that in the model ( ^.15|) the operators 
q^ = ql^~~^^ and qk = qk'~^^ are related as ql = qk for arbitrary k. This identity together with 
leads to a similar relation between possible values of the quantum numbers 



Qk = Qk 



{k 



and as a result the energy of the compound Reggeon states is equal to 



E„ 



{q)+en{q*)) 



27r 



ReSniq) 



(5.24) 



(5.25) 



in accordance with (|5.21|) . We can find possible values of q2 and q2 in (|5.24| ) using the inter- 
pretation ( [4.17|) of the corresponding operators as the Casimir operators of the SL{2,€) group. 
Relation (|5.4|) implies that the eigenvalues q2 and q2 are the same for the Reggeon hamiltonian 



and for the hamiltonian (|5.15|) . In the model ( ^.15|) , the operators and act on the 



total quantum space t 



0,2 , 



I which can be decomposed into the direct sum of irreducible 
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principal series representations, (Bu,m.t'^''^"^, of the 5*17(2, C). Then, for each t'^'^™ the Casimir 
operators q2 and q2 have the following eigenvalues 



q2 = -h{h - 1) , q2 = -h{h - 1) 



1 + m 



1 — m 



— IP , 



— IV 



(5.26) 



with integer m and real v. One checks that these expressions are in accordance with ( |5.24|) . 

Invariance of ( |5.22D and ( |5.23|) under the conformal transformations ( p.4|) implies the fol- 
lowing transformation properties of the Reggeon wave function. 



Xn,{a}{{zu Zi}; Zq, Zq) Xn,{a}iWv ^'i}'^ 4' 4) = (^^^O + d) {cZq + d) Xn,{c}{{Zi, Zi}; Zq, Zq) 

with conformal weights h and h defined in (|5.26| ). This relation can be rewritten in terms of 
generators of the conformal group as follows 



Y.^k + do] Xn. 



{a} 



Xn,{a} Xn,{a 







\k=l 



\k=l 



Vfc=l 



and analogously for the antiholomorphic generators. Here, the first equation means that 
Xn,{a}{{zk, Zk}; Zq, Zq) depends on the distance Zko = Zk — zq to the center of mass, and the 
two remaining equations can be represented as 



Xn,{a}{{Zi, Zi}; 0, 0) = -hXn,{a}{{Zi, Zi}; 0, 0) 



Xn,{a}i{z^,z,};0,0) = (5.27) 



with S"" the total spin for s = 0. One can show that the wave functions of the Reggeon 
compound states, ( p.lO|) and ( |5.13| ), satisfy (|5.27|) . 

Using relations ( ^^221) and ( ^^231) we can give a meaning to the decomposition of the transi- 
tion amplitude ( |2.12| ) over the Reggeon states. Comparing ( ^.23| ) with analogous decomposition, 
I = E{a} \x){x\: used in (|2.12|) , we obtain. 



IX) = Xn,{q}{{Zk, Zk}; Zq, Zo) 



{X\ = {Xn,{q}{{Zk,Zk};Zo,Zo)) {Zi2Zi2...ZnlZnl) ^ 



Let us consider now the state {x\Tn'^ which enters into ( ^.12| ) and ( ^.14| ). We recall that 
operator T^^ describes the free propagation of n Reggeons, eq. ( |2.9| ), and therefore the state 
(x|T^°^ can be represented as 

(Xirf) = {x\{PlPl-PnPny' = {{qnqnr'XnM} 

with operators g„ defined in ( f4.16| ). Since Xn,{q} diagonalize g„ and the operators can be 
replaced by their corresponding eigenvalues, g„ and g*, respectively. After substitution of the 
last two relations into ( |2.12|) we find transition operator as 



Tn{{z, z}, {z, z'}; uj) = Y^ I dzodzo 



Xn,{q}{{z, Z}; Zq, Zq) {xn,{q}{{z' , z'}; Zq, Zq) 



uo-E. 



■n,{q} 



. (5.28) 



Here, the sum goes over all possible compound n— Reggeon states propagating in the t— channel. 
Integration is performed over the position of the center of mass bo = {zq, zq) of the states and 
after Fourier transformation / d^boexp{ibo ■ qo) it can be replaced by integration over the total 
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transverse momentum go carried by the Reggeons in the t— channel Momentum conservation 
fixes go to be the transferred momentum g in the scalar products entering (|2.4| ) and ( p.l4| ). As 
a result, for the residue factors ( |2.14| ) we obtain the following expressions 



PaHq^) = fdzodzo e^'^-'^fdzdz $a({z, z})xn,M({^, z}; Zo, Zo) , Pi'\q') 



(5.29) 

We conclude from ( |5.21|) and ( |5.29|) that, in accordance with our expectations, the energy of 
the Reggeon states and the product of the residue factors are real in ( [2.131 ). 



6. Solving the Baxter equation 



To find the explicit expressions for the spectrum of the Reggeon hamiltonian, ( [5.25| ) and ( |5.10| ), 



we have to solve the Baxter equation ( |5.9|) for an arbitrary n and for fixed set of the quantum 
numbers {g^}. The next step will be the identification of the values of h and {qk} corresponding 
to maximum eigenvalue -E™^^, which determines the dominant contribution to the hadronic 
scattering amplitudes ( |2.13| ) and ( P^.15| ). To understand better this procedure we first consider 
the simplest case n = 2. 



6.1. Special case: 



n 



For n = 2 the Reggeon hamiltonian (|3.9|) is equal to 2{Hi2+Hi2) with the two-particle XXX spin 
s = hamiltonians H12 and H12 given by ( [4.5|) . They depend on the operators J12 and J12, which 
are related to the Casimir operators of the conformal group as Ji2(^i2 + 1) = —^2 = — -2i2"9ic?2 
and Ji2(</i2 + l) = —q2 = —zf2did2- Possible eigenvalues of these operators were found in (|5.26| ) 
and the corresponding eigenstates define the wave function of the compound state of n = 2 



Reggeons, the BFKL Pomeron [p7? 



Zo, Zo) 



l + m 
2 



. / zu A~"'V 



Zl2 



V ^10^20 



(6.1) 



where the constant N^^u = '^{^'^ + ^^1^/4) fixes normalization of the state. To find the energy 
corresponding to this state one can use either the definition ( [4.5|) , or apply the relations (|3.11|) , 
(|5:T6| ), (|5l^ ) and (^^, 



E2{m, v) = ^^2 Re 



TT 



/ 1 + Iml ^ 



(6.2) 



Thus, the energy of n = 2 Reggeon state is finite and real. The maximum value of the energy 



^2(0,0) = ^4 log 2 



TT 



(6.3) 



corresponds to m = = 0, or equivalently to h = h = ^. One notice that the expression (|6^ ) is 
invariant under replacement h —>■ h* and h ^ h, and the maximum value of E2 corresponds to 
the values of h and h which are not changed under this transformation. This simple observation 
will be generalized in Sect. 6.3 for the spectrum of the Reggeon hamiltonian with n > 2. 
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6.2. Properties of the Baxter equation 

After one introduces the notation (|5.1CI| ) for the function Qn{^) = A"(5n(A), the Baxter equation 
(|5.9| ) takes a form of the discrete one-dimensional Schrodinger equation with a singular potential 

Qn(A + + Q.(A - ^) - 2Q„(A) = (^-^^^^^ + ? + •■• + I^) ^"^^^ (^-^^ 

where we replaced q2 by the expression (|5.26|) in terms of the conformal weight h. We notice 
that ( |6.4| ) is invariant under the replacement h ^ 1 — h and we may restrict the values of the 
conformal weights to be 

Re/i>l/2 (6.5) 

or equivalently m > in ( p. 26 ). Taking a naive continuum limit in ( |6.4|) , A — > oo, we find that 
for finite h and {q} the equation has two independent solutions corresponding to the different 
asymptotic behavior of the solutions 

A"g„,(A) CiA'' + C2X^-^ (6.6) 

with Ci and C2 arbitrary constants. One of the solutions grows at infinity while the second 
one vanishes. To make a choice between them, we compare ( |6.6|) with the expression (|5.14| ) 
for Qn{X) given by the algebraic Bethe Ansatz. We find that C2 = for integer conformal 
weights h > n. This suggests to impose the following asymptotic behavior on the solutions of 
the Baxter equation (^.9]) ,["] 

Q„(A) A''-" (6.7) 



with the conformal weight h in the fundamental domain (|6.5|) . 

Suppose we know the solution of the Baxter equation, Q„ = (5n(A; /i, {g}), which obeys 
(|6.7| ). Then, replacing /i— >1 — /lorA— ^ —A and qt —>■ {—)'^qk in (|6.4|) we find that equation 
takes the original form, which means that up to inessential factor 

g„(A; h, {qk}) = g„(A; 1 - h, {q^}) = g„(-A; h, {{-fqk}) ■ (6.8) 

We substitute this identity into expressions for the energy and for integrals of motion, (|5.10|) , 
and find that 

en{h, {qu}) = £„(! - K {qk}) = £nih, {H^k}) (6.9) 

and 

Mh, {Qk}) = Xfc(l - h, {qk}) = i-)'lk{h, {H^qk}) • 

The last part of these relations can be understood from the properties of the model under 
reparameterization of the Reggeon coordinates, Zk — > Zn-k+i, with k = l...n. We check that 
the Reggeon hamiltonian (|3.9|) is invariant under this transformation while the operators qk, 



defined in ( |4.16| ), change sign for odd fc. 

Let us consider the special values of the quantum numbers, qk = qk+i = ■•• = Q'n = with 
n the number of Reggeons and A; > 3. Then, the solution of the Baxter equation (|5.9| ) for n 
Reggeons can be expressed as 



g„(A; h, gg, gfc-i, 0, 0) = A-"+^-^gfc_i(A; h, gs, Qk-i) (6.10) 



^°This condition fixes the ambiguity of ( |6.4| ) those solutions can be defined up to factors exp(±27rAfc) for 
k = which lead to the exponential growth of at infinity. 
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where Qk-i{X; h, q^, qk-i) is the solution of the Baxter equation (|5.9| ) for /c — 1 Reggeons 
with the same quantum numbers h and gs, ... , qk-i- To understand this identity we substitute 
it into ( [5.10| ) and find the relation between the eigenvalues of the holomorphic hamiltonians for 



different number of Reggeons 

en{h, gg, qk-i, 0, 0) = Ek-iih, q^, g^-i) • (6.11) 

As it follows from the definition, (|4.16|) , the operators are given by a sum of differential 
operators and each of them acts only on k holomorphic variables from the total set of n Reggeon 
coordinates zi,...,Zn- Therefore, to satisfy qk = qk+i = ■■■ = qn = one can choose the 
eigenstate to be a function of only k — 1 holomorphic Reggeon coordinates, that is to be the 
holomorphic part of the compound state of k — 1 Reggeons. Using ( ^.28|) we find that these 



states give infinite contributions to the transition operator and therefore should be excluded 
from the spectrum of n— Reggeon hamiltonian. This means, that solving the Baxter equation, 
(|5.9| ) and (|6^ ), we should satisfy condition 

Qn^O. (6.12) 



In particular, for n = 2 this condition together with (|5.26|) and ( |6.5|) implies that h ^ 1 and 



6.3. Maximum eigenvalue of the Reggeon hamiltonian 



Let us consider expression ( ^.251) for the energy of the n— Reggeon compound state as a function 
of possible sets of the quantum numbers En = En{q)- Then, it follows from ( |5.25| ) and (|6.9|) , 
that this function has the following properties 

Eniq) = En{q*) = Er,{-q) , 

where q = (g2, gs, gn) and — g = (g2, — gs, (— )"gn) and where complex conjugation acts 
in g* on all elements of the set. This means that the spectrum of the Reggeon hamiltonian is 
degenerate with respect to quantum numbers {g^}. Hence, if we assume from the very beginning 
that the eigenstate of the Reggeon hamiltonian corresponding to maximum eigenvalue is not 
degenerate, then the values of the corresponding quantum numbers can be fixed as|^ 



q2k 



(l2k 



q2k 



(l2k 



q2k+l 



q2k+l 



0, 



In particular, for quantum numbers g2 and g2 defined in (|5.26|) this relation leads to 

1 



92 



92 



for h = h = 1/2 



4 

max t: 

which indeed corresponds to the maximal energy ( |6.3D in the special case n = 2. 



(6.13) 



(6.14) 



Changing the sign of the Reggeon hamihonian we can reformulate this condition as uniqueness of the ground 
state of the model. 
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6.4. Ansatz for the Baxter equation 



To solve the Baxter equation ( |5.9|) we choose the following ansatz [|T9| 



(6.15) 



where the function Qn{z) is integrated along the closed path Cz in the complex plane. 
Integrating by parts k times we get the equivalent representation 



Qn{X) = X~'j ^{Z 



d \ 

-iz{l - Z)—] Qn{z) 



(6.16) 



which allows us to replace the original functional equation ( p.9| ) by the n— th order differential 
equation on the function Qn{z) 



z{l-z) 



dz 



n-2 



+ z{l-z)Y.t''-''q^Jz{l-z) 



k=0 



dz 



Quiz) = 



(6.17) 



with q2 = —h{h — 1). It is important to notice that the integration by parts in ( p.l6| ) leads to 
the integrals like J^^dz^^...) where (...) denotes the product of Qniz) and z~'^^~^{z — 1)*'^^-'^ 
or their derivatives with respect to z. The choice of the integration contour is fixed by the 
condition that these additional terms can be neglected, or equivalently that after encircling the 
path Cz the function (...) comes to the initial point with the same value. To this end we notice 
that z~*'*'~^(2; — 1)*'^"^ has a nontrivial monodromy at z = and z = 1 and the same points are 
the singular points for the differential equation (|6.171) . 

The asymptotics ( |6.7| ) of the solution of the Baxter equation corresponds to the behavior 
of the function Qn{z) at infinity 

Q^^.y^ .l^^n+l ^ (6.18) 

which may be considered as a possible additional condition on the solutions of the differential 



equation (|6.17|) 



6.5. Solution for n = 2 

For n = 2 the differential equation (|6.17| ) has a form of the hypergeometric equation 



dz 



z{l-z)-^ + h{h-l)\Q2{z) = 



(6.19) 



and its general solution is known as ^6 



Q2{z) = const. / ^ w^'-'iw - l)''-\w - z) 
JCu, 2m 



-h 



(6.20) 



where the integration path Cw is fixed by the condition that the integrand should be a single 
valued function on C^, that is after encirchng around it should have the same value. We 
substitute (|6.20|) into ( |6.15D and find the general solution of the Baxter equation for = 2 as 



(52(A) = const. / -f^ z~ 

JCz 2m 



'^~Hz-ir^-^ 



Cn, 2m 



w — z) 



(6.21) 
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This expression depends on the integration paths Cz and Cw Their choice depends on the 
values of A and h. To show this let us consider the special case of integer positive conformal 
weight h. For h = Z+, the integrand in ( |6.21| ) has a pole at w = z and for the integral over 
w do not vanish the contour in ( |6.2(]| ) and (|6.21 ) should encircle the point w = z. The 
result of taking residue aX w = z leads to the well-known Rodrigues's formula for the Legendre 
polynomials^ 

Q2{z) = Ph-i{1 - 2z) = —l—^{z{l - z)f-' = F{h, l-h;l;z). 



{h - 1)\ dz''-^ 



The integral over z in (|6.21|) has two branch points at 2; = and z = 1. For the integrand 
to be single valued, the contour should encircle both branch points in the same direction. 
Then, we deform the integration path to be close to the segment [0, 1] and get the solution of 
the Baxter equation for integer h as 

■h+i 



Q2{X;h) 



sinh(7rA) 



dz z-'^-\i-zy^-'Ph-i{i 



2z) 



(6.22) 



where the constant has been choosen for the solution to be real for h = 2+. The Legendre 
polynomials Pk{x) for k = 0,1,... form a system of orthogonal polynomials on the interval 
— 1 < X < 1 with the weight 1. They have the following properties p7[ 



and 



.(l-.)^P.(l 



Pk{x) = P„fc_i(x) , Pk{-x) 

d ^ _ k{k + l) 



-lYPkix) 



2z) 



Using (|6.23| ) we verify that 
g2(A;/i) = ^'-'^Q2(A;l-/i) 



2(2A;+ 1) 
Q2{-\\h)-. 



fc+i 



1 - 2z) 



Pi 



k-l\ 



2z) 



(6.23) 

(6.24) 



Q;{X;h) = Q,{X;h). 

(6.25) 

Applying identity ( |6.24| ) together with (|6.16|) we find that Q2{X; h) satisfies the recurrence 
relation 

[g2(A;/i + l) + g2(A;/i-l) . 



Ag2(A;/i) 



(6.26) 



2{2h - 1) 

Using expansion of the Legendre polynomial Ph-i{2z — 1) in powers of z we obtain the repre- 



sentation for the solution of the Baxter equation for integer conformal weight /i as 



Q2{\\h) 



fc=i 



k T{h + k)T{k-i\) 

' {k\ fr{h-k) r(i -zA) 



z^/i(l - h) 3F2 



'l + h,2-h,l 
2,2 



iX 



For lowest values of the conformal weight h the explicit expressions for g2(A; h) are 



(6.27) 



Q2(A;2) 



g2(A;3) = -6A, g2(A;4) = -2 + 10A^ g2(A;5) 
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A 



35 



A^ 



(6.28) 

As it was expected from the algebraic Bethe Ansatz, ( |5.14|) , they are polynomials of degree 
h — 2 in X and one checks that the roots of g2(A; h) do satisfy the Bethe equation 



"'^^We should notice that there is a second solution of ( ^.19| ), kn own as the Legendre function of the second 
kind 1 36, 3^. However, it does not s atisf y the additional condition ( 3.1g ) and leads to the solution of the Baxter 



equation with asymptotic behavior (6.6) for Ci = 0. 
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6.5.1. Relation to orthogonal polynomials 

The remarkable property of the expression (|6.27| ) is that being taken in the form of 3F2 hy- 
pergeometric infinite series for arbitrary complex h and A, the expression for Q2(A; h) becomes 
a polynomial of finite degree either for h = m or for X = im with m integer. This suggests 
that there is a deep relation between solutions of the Baxter equation for n = 2 (and higher 
n) and some systems of orthogonal polynomials. The latter relation was explicitly specified by 
V. Spiridonov who has observed that solutions of the Baxter equation for n = 2 are tied to the 
Hahn class of orthogonal polynomials. 

Indeed, if one sets A = —i{m + 1) with m = 0, 1, ... and h = ix + 1/2, then the Baxter 
equation ( |5.9| ) takes the form of three term recurrence relation which defines the continuous 



dual Hahn polynomials of the argument (see, e.g. ||38||), 

i,a + ix,a — ix 



■I/O f \ / 1 \ / \ ( TTX/ .J Gj I IjOb .J Gj I/Ob \ . _ . 

m(x ;a,6, c) = (a + 6)m(a + c)m 3F2 ;1 , a = 3/2, 6 = c = 1/2 . 

y a + b,a + c J 

These polynomials are orthogonal on the interval < x < 00 with respect to a positive measure 
and they coincide up to an inessential factor with the solution Q2{— 1(171 + 1); ix + 1/2). 

Due to the well known property of duality of classical orthogonal polynomials the same 
solution of the Baxter equation determines different system of orthogonal polynomials for which 
A becomes continuous variable while h plays a role of integer index. Indeed, putting h = 171 + 2 
with m = 0, 1, ... and A = a; in ( |6.26| ) we find that ( |6.26| ) takes the form of three term recurrence 
relation for continuous symmetric Hahn polynomials ||38|| , 



, n .m (-m,m + 2a + 2b-l,a-ix \ 

which are orthogonal on the interval —00 < x < 00 and which are proportional to Q2{x; m + 2). 

Identification of Q2(A; m + 2) as Hahn polynomials implies, first, that the solutions, {A} = 
(Ai,... ,Xk-2), of the Bethe equation ( |5.7|) and (|5.14| ) have the properties of the roots of or- 
thogonal polynomials [^. Namely, all roots {A} are real and simple, between two consecutive 



zeros of Q2{X; k) there is exactly one zero of Q2(A; k + 1) and at least one zero of Q2{X; k + m) 
for each m > 0. 

The relation between orthogonal polynomials and solutions of the Baxter equation is a 
general property of the Baxter equation for XXX magnet of an arbitrary spin, including the 
case of compact SU{2) spins. In particular, some solutions of the Baxter equation ( |5.9| ) for 
n = 3 can be identified with the Wilson polynomials. The detailed consideration of these and 
other analytical solutions of the Baxter equation for generalized Heisenberg magnet models on 
the basis of methods developed in |39] will be given elsewhere . 



6.5.2. Eigenvalues for n = 2 



Let us use the solution ( |6.27| ) to obtain the spectrum of the Reggeon hamiltonian for n = 2 



which has been already defined in sect. 6.1. We start with the expression for the energy, ( |5.10D 



and evaluate Q2(A; h) and its derivative at A = — i while the same quantities at A = i can be 
found from (|6.25|) . We use ( |6.27| ) to get for h = 



Q.,-. . /|(-l)^^£|;i±|£(i=^ . ^,(1 - ft) (6.29) 
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since only one term with k = 1 survives in total sum. For the derivative we obtain 



-^^<52(-A-z;/i) 



h-l 



T{h + k) 



u=2 (k - mh - k)v\k + 1) • 

To evaluate this sum we rewrite it as a contour integral in the complex j— plane 



A=0 



TT 



dj 

2m sin(7rj)(j-l)r(/i-j)r2(j + i: 



where integration contour encircles half of the real axis right to j = 2. The integral is defined 
by singularities of the integrand in the complex j— plane. A simple analysis shows that the 
integrand has poles at j = 0, 1, 2, ...,h — l while singularities at j = —1, —2, . . . are compensated 
by r^(j + 1) in the denumerator. Since the points j = 2,3, h — 1 are inside the integration 
contour, the integral is given by the residue at j = and j = I, 



-2&h{l-h) [i){h)-i){l) - 1] 



(6.30) 



A=0 



with h positive integer conformal weight. Now, we substitute ( |6.27| ) into ( |5.1CI| ), take into 
account ( |6.29|) , ( |6.30| ) together with (|6.25|) and obtain the following expression for the energy 
of the holomorphic Reggeon hamiltonian 



e2{h) = -A[ij{h)-i,{l)] 



for /i = Z+ . 



(6.31) 



Expression ( |6.31 ) is valid for positive integer conformal weight h and in order to find the energy 
for negative integer h we have to use the symmetry ( |6.9| ). In particular, the conformal weight h 
corresponding to the antiholomorphic hamiltonian can be found from (|5.26|) ash = l — h* and 
for positive h it becomes negative or zero. In general, for an arbitrary complex h the energy of 
the n = 2 Reggeon state is given by ( |5.25| ). We replace e2{h) in ( |5.25| ) by its expression ( |6.31 



for integer h and analytically continue the result to complex values of the conformal weight. 



O) 



This expression coincides with the energy 



of the Reggeon hamiltonian at n = 2. 



6.5.3. Eigenstates for n = 2 



For n = 2 we obtain from (|5.13| ) the expression for the holomorphic wave function of n = 2 
Reggeon state as 

iPh{zi, Z2, Zq) = -z{^ {iS'Y^'^ Q2{^\\ h) „ „ , (6.32) 



2 2 ' 
^10^20 



where the operator xi was defined in ( |5.12|) and S = —{di + 82). Substituting ( |6.28|) into 
(|6.32|), we find that 



^2 



<— ) 

V ^10^20/ 



\z10z20y 



\Z10Z20y \ ^10^20/ 
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To find (fh for an arbitrary h we introduce the notation for the holomorphic part of the wave 
function (|5.16|) 

HS-)''-^ Q2{xi;h) -J-^ (6.33) 



"^12 Vh 



^10^20 



and use the property ( |6.26 ) of Q2{xi; h) together with (|6.32 ) to transform it into the recurrence 
relation 



did2Zi2(t)h 



h{h-l] 



2{2h - 1, 

One can check that the solution of this equation is 



+1 - {di + 82 



h-2 
Z\2 



(^10^20) 



where = hNh-i- Substituting (ph into (|6.33| ) we obtain the holomorphic wave function of 
the n = 2 Reggeon state in the form ( |6.1| ). Thus, for n = 2 the solution of the Baxter equation 
(|6.27| ) allows us to reproduce the spectrum of the Reggeon hamiltonian. 



6.6. Solution for n > 3 



Analyzing the differential equation ( 6.17|) for n > 3, we follow the same procedure as for n = 2. 
Namely, we start with integer positive conformal weight h and try to find solution of (|6.17|) and 
(|6.15| ) which can be analytically continued for complex values of h given in ( ^.26|) . However, 
for n > 3 the solution of (|6.17|) depends on additional quantum numbers q^, which makes 

analysis more complicated. 

Let us introduce notation for the derivative 



Qn{z) 



-iz{l — z) 



dz 



n-2 



Qn{z) 



Applying operator \ z{l — z)j^ 
equation for Qn{z) as 



n-2 



".{l — z)) ^ to the both sides of (|6.17|) we find the differential 



zil-z) 



dz 



n-2 



n-2 



^z{l-z)^ + Y.^^-'<ln- 



zil-z) 



dz 



Qn{z)=0 



(6.34) 



Comparison of this equation with (|6.19|) and the property ( |6.24|) of the Legendre polynomials 
suggest us to look for the solution of ( |6.34D as a linear combination of the Legendre polynomials 



(6.35) 



k=l 



where summation is performed over positive integer k and the coefficients Ck = Ck{h, {q}) have 
to be found from ( |6.34| ). Once we know these coefficients, we use ( |6.15| ) and ( |6.17] ) to find the 
solution of the Baxter equation as 



Q.(A; h, {q}) = A^-" ^2(A; k) Ck{h, {q}) 

k=l 



(6.36) 
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with Q2{X',k) the solution of the Baxter equation for = 2 defined in ( |6.27| ). Substituting 
this expansion into (|5.10| ) , we find expression for the energy of the holomorphic hamiltonian in 
terms of the coefficients Ck 



Snih, {q}) = -4 ^ - V^(l)] Ckih, {q}) 



k=l 



(6.37) 



with normahzed coefficients Ck defined as 



a = Re 



i''k{k - l)Ck 



'Em>i i'^rn{m - 1)C„ 



(6.38) 



fc=i 



For n = 2 we have = 5k,h and only one term with k = h survives in the sum ( |6.36| ) and 
(|6.37| ). For n > 3 we substitute (|6.34|) into (|6.17|) and equate to zero coefficients in front 
of Legendre polynomials. To this end, we notice from ( |6.24| ) that acting on the Legendre 
polynomial, Pfc(l ~ 2^;), the operator z{l — z)-j^ shifts its index, while operator -j^z{l — z)-^ is 
diagonal. Hence, in the basis of the Legendre polynomials, \k) = Pfc_i(l — 2z), both operators 
can be represented as infinite dimensional matrices 



m,k 



■k~m—l 



{m\z{l 



k{k 



2{2k - 1 



m,fc+l 



m,k—l ) ) 



m,k 



'{m\—z{l — z)-^ + h{h — l)\m) 



{k + h - l){k - h)6m,k 



(6.39) 



In terms of these operators, the resulting equation for the coefficients Ck has the following 
matrix form 

W{q) = [u^-^V + j^qk f/"-'^) , W{q)\C) = , (6.40) 

\ A:=3 / 

where \C) denotes infinite dimensional vector with components Ck and the term with q2 has 
been included into the definition ( |6.39|) of the operator V . 

The condition ( 6. IS ) leads to the asymptotic behavior of the function Qn{z) 



Qr, 



for h = ^. 



(6.41) 



Let us look for the solutions of ( |6.34| ) and ( |6.35| ) which are finite sums of the Legendre polyno- 
mials. Then, from the asymptotics of the Legendre polynomials at infinity, Pn{z) ~ -z""^, and 
from comparison of ( |6.41| ) and ( |6.35D we find that 



c 



h+1 



c 



h+2 



(6.42) 



and nonzero components of \C) form a /;,— dimensional vector. It is clear from ( 6.401) that the 
resulting finite-dimensional vector is the eigenstate of the hxh minor of the infinite-dimensional 
matrix ( |6.40| ) corresponding to the zero eigenvalue. This implies that the determinant of the 
minor should vanish which leads to the condition on the quantum numbers qk entering into 
( |6.40| ). Moreover, as we will show below, there are n — 3 additional conditions on qk which 
follow from the relations ( |6.42 ) and hence the total number of the independent constraints 
turns out to be equal to the number of gfc's. This means, that for the Baxter equation ( |5.9| ) to 
have polynomial solutions, the quantum numbers qk have to be quantized. 
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6.6.1. Special case: n = 3 

Let us consider tlie explicit form of (|6.40|) for n = 3, 

fk+iCk+iils) + fk-iCk-iiQs) 



qsCkiqs 



where the coefficients fk depend on the conformal weight h and are given by 

k{k-l) 



fk 



2{2k - 1) 



{k-h){k + h~l) 



fi = fh = 0. 



(6.43) 



(6.44) 



The last identity ensures the trancation of the recurrence relations ( |6.43D at k = h — 1. Intro- 
ducing the notation for Wh to he h x h minor of the matrix W defined in ( |6.39|) and for \C)h 
to be /i— dimensional component of the vector \C), 



( qs -f2 

-/l <?3 -/s 

-/2 qs 



\ 







\ 



— fh 

-fh-1 gs J 



( Ci \ 
[a] 



(6.45) 



we rewrite original matrix equation ( |6.4CI| ) as 

detiy;,(g3) = 0, 



(6.46) 



Here, the first equation gives us possible values of the quantum numbers gs for positive integer 
conformal weight h, while the second one determines the coefficients entering into solution of 
the Baxter equation ( |6.36| ) and energy of the Reggeon hamiltonian, ( |6.37| ) and ( |6.38| ). It can 
be easily shown from ( |6.45|) that 



detH^,(-g3) = (-l)'detl^,(g3). 



(6.47) 



For even h this equation has h — 2 nontrivial solutions, q^ ^ 0, and a degenerate solution = 0, 
which being substituted into ( |6.43| ) leads to Ck = 0. For odd h it has h — 3 solutions for gs 7^ 
(see fig. ID and a solution gs = corresponding to the degenerate states discussed in Sect. 6.2. 
Hence, for fixed h we have 2[/i/2 — 1] nontrivial solutions for which satisfy (|6.12| ) and which 
determine the solutions ( |6.36| ) of the Baxter equation, or equivalently 2[h/2 — 1] eigenstates of 
the holomorphic Reggeon hamiltonian. Although it is not difficult to solve (|6.43| ) or (|6.46| ) for 
lowest values of /i = Z+, the general analytical formula is not available yet. 

Since detWh^qs) is a polynomial of degree h in q^ with real coefficients and certain parity 
( |6.47| ), its zeros may be either real or pure imaginary. It turns out that only the first possibility 
is realized 

Im = , for /i = Z+ . (6.48) 

To understand this property we notice that the coefficients Ck^qs) satisfy the three-term recur- 
rence relations ( |6.43| ), which define a system of generalized orthogonal polynomials in continuous 
variable qs-^ Condition Ch+iiqs) = is equivalent to ( |6.46|) and it implies that the possible 



^■^The recurrence relation similar to ( 3.43 ) and the interpretation of its solutions in terms of orthogonal 
polynomials have been also proposed in [ ll[. 
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values of gs, being roots of orthogonal polynomial, are real and simple. We conclude that each 
term in the sum (|6.36|) which determines the solution of the Baxter equation for n = 3 and 
/i = Z has a form of the product of two orthogonal polynomials, Q2{^', k) in continuous variable 
A and Ck{qz) in continuous variable ^3. 

One can find Ck from ( |6.43| ) and ( |6.46| ) up to some constant which can be fixed by putting 
Ci = 1. Then, it follows from ( |03D and (^M ) that Ck are real and Cki-qs) = i-)'"^^Ck{q3)- 



As a consequence the solution of the Baxter equation ( |6.36| ) and the energy of the holomorphic 
hamiltonian (|6.37| ) are also real for integer conformal weight h and real A. 

Here, we give few examples of the solutions of the Baxter equation, Qsi^X; h, q^), for lowest 
values of the conformal weight h and positive q^, 

Qs{X;A,2V3) = -6 - 6 VSA , 

Qs{X; 5, 6^3) = -f - lOVSA - f A^ , 

A - 30A^ 



Q3(A;6,4V30) 
g3(A;6,4v^) -- 



5V30 



15 
2 



35^3 
4 



A+f A2+ 25^A^ 



QsiX; 7, 6V23 ± 6v^) = -^i^^ - 2h/2M^X - m±I^X^ - 2lV23|119V3^3 _ 63±49v^A 



For h < A the Baxter equation has only degenerate solutions ( |6.10 ) corresponding to qs = 0. 
As a nontrivial check of these expressions we verify that, similar to the case n = 2, the roots 
of Q3(A;/i, gs) satisfy the Bethe equation ( p^) for n = 3. Moreover, we find that the roots 
{Afc} are real and simple. This should be compared with the situation for spin s = 1/2 XXX 
magnet, where the solutions of the Bethe equation ( p.7[ ) form strings in the complex A— plane 
symmetric with respect to real axis [|^, |2^, |2^, |2^ . 

Unfortunately, we do not have analytical expressions for Q^IX; h,q3) similar to that for 
n = 2, eq. (|6.27|) . For lowest values of h the results of our calculations of the solution of the 
Baxter equation and the energy of the holomorphic Reggeon hamiltonian are shown in fig. 
For fixed /i = Z_|_ the energy e^^h, q^) turns out to be a decreasing function of ^3 and for — > 
it approaches the value e2{h) = —4:[ip{h) — ip{l)] in accordance with ( |6.11D . 



6.6.2. Generalizations for n > 4 



Let us consider the properties of the equations (|6.40|) and (|6.42|) for n > 4. Using the explicit 
form of the matrices ( |6.39D one finds that for an arbitrary n the coefficients Ck satisfy the {2n — 
3)— term recurrence relations similar to ( |6.43|) which relate the coefficients Ck-n+2, ■■■,Ck+n-2 
and involve the quantum numbers q3,...,qn. If we again define Wh as a /i— dimensional minor 
of the matrix ( |6.40| ), then the relations ( |6.46| ) are still valid. The first equation in ( |6.46D gives 
the relation between g3,...,g„, while the second one allows us to find the coefficients Ck as 
functions of these quantum numbers. We notice that for = 3 it is enough to satisfy Ch+i = 
in order to get Ch+2 = Ch+3 = ... = from the three-term recurrence relations ( |6.43|) . For 
n > 4 the recurrence relations for Ck involve {2n — 3) terms and as a result equation ( |6.42| ) 
becomes equivalent to Ch+i = ■■■ = Ch+n-2 = 0. These relations together with ( |6.40| ) lead to 
the following conditions 



W^\q)Ck = Q, (6.49) 
h + n — 3. Once we know expressions for \C)h from ( |6.46| ), 



with /c = 1, h and m = /i + 1, 
we substitute them into ( |6.49|) to get {n — 3)— additional conditions on q^ 



, qn which together 



33 



with (|6.46|) completely fix allowed values of these quantum numbers. The problem becomes 



pure algebraic and we summarize the properties of the solutions in figs.^|TO|. 

It turns out that for ri > 4 and integer positive conformal weight, h = Z_|_, the Baxter equa- 
tion (|5.9|) has nontrivial polynomial solutions only for h > n since for h < n one can not satisfy 
condition ( |6.12| ). Similar to ( |6.48| ), the quantum numbers q3,...,qn which one finds solving 



(|6.46| ) and (|6.49|) are real and simple. Then, the solution of the Baxter equation, Qn(A; h, {q}), 



becomes a polynomial of degree h — n in spectral parameter A with real coefficients. The roots 



Ai, Xh-n of Qn{X; h, {q}) satisfy the Bethe equation (|5.7|) . We find that for each solution of 
the Baxter equation, the roots are simple and real 

ImAfc = (6.50) 

with k = 1,2, h — n. Being taken together, the roots corresponding to all solutions of the 
Baxter equation for fixed n belong to the interval of real axis in the complex A— plane symmetric 
with respect to 0. The length of the interval increases as n grows. 

The property ( |6.50| ) allows us to find lower and upper bounds on the energy of the Reggeon 
hamiltonian, £„(/i, q), corresponding to the polynomial solutions of the Baxter equation. Using 
( 6.50 ) we find from ( ^.8|) the following relation 



2h<en{h,q)<-2n. (6.51) 



One checks, that for n = 2,3 and 4 this estimate is in agreement with our results, ( |6.31| ) and 



fig. p|-p!0|. Moreover, one can find using ( |5.10|) and ( ^.14|) that the maximum value of energy 



Sn = —2n, is achieved for the trivial solution of the Baxter equation (|5T 
Qn{X) = const. , h = n, q2k = — 7^ , q2k+i = , Lpn 

'nO 



{2k)\{n-2kyr ^^^-^ ^' z.Vio-^. 



which is in agreement with ( |6.13| ). Notice that for odd n this solution is degenerate since it 
does not satisfy (|6.12|) . 

As we have seen in Sect. 2, the Regge asymptotics of the scattering amplitude is dominated 
by the n— Reggeon states with maximum value of the energy ( 5.25| ). The relation (|6.51|) implies 



that using the polynomial solutions of the Baxter equation we cannot exceed the bound for the 
energy En = —^^n. However this value cannot be maximal because for n = 2 the maximal 
value of the energy, ( |6.3|) , is positive. This is in accordance with our expectations, (|6.13| ) and 



( |6.14| ), that the maximal value of energy corresponds to the conformal weight h = 1/2, while 
polynomial solutions exist only for integer h > n. Nevertheless, the reason why we might be 
interested in solving the Baxter equation for integer conformal weight h is that, similar to the 
situation for n = 2 described in Sect. 6.5, we hope to find expressions for the spectrum of 
the Reggeon hamiltonian which can be analytically continued for arbitrary h. The analytical 
expressions for n = 2 and the properties of the Baxter equation for n > 3 (see figs.^ p!0| ) 
indicate that the polynomial solutions of the Baxter equations should be related to new systems 
of orthogonal polynomials . 



7. Conclusions 

In this paper we studied the asymptotic behavior of hadronic scattering amplitudes in the 
Regge limit (|1.1| ) of high energies and fixed transferred momentum. The consideration was 
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performed in the generalized leading logarithmic approximation which preserves unitarity of 
the S*— matrix and which requires the resummation of infinite series of nonleading logarithmic 
corrections, log"^ s, to the scattering amplitude. Although this approximation is rather crude, 
e.g., it does not allow us to fix the argument of the coupling constant in ( p.l3| ) and (|2.15|) , it 
gives the first order approximation to the exact result. 

After resummation to all orders of perturbation theory, the scattering amplitude has a 
Regge behavior which can be associated with the singularities of the partial waves in the com- 
plex angular momentum plane. In the generalized leading logarithmic approximation, due to 
gluon reggeization, longitudinal gluonic degrees of freedom become semiclassical and nontrivial 
interaction occurs only in two-dimensional space of transverse gluon momenta. As a result of 
this interaction, reggeized gluons form color singlet compound states - Pomeron, Odderon and 
higher Reggeon states, which satisfy the Schrodinger equation with the Reggeon hamiltonian 
depending on the color charges of gluons and their transverse momenta. The Reggeon hamilto- 
nian describes the pair-wise interaction between reggeized gluons and it obeys the remarkable 
properties of holomorphic separability and conformal invariance in the impact parameter space. 
The properties of the compound Reggeon states govern the large— s behavior of the total cross 
section, (|1.3| ) and ( p.l3|) , and small— x asymptotics of the structure function of deep inelastic 
scattering, ( p. 151 ). 

In the impact parameter space, holomorphic and antiholomorphic gluonic degrees of freedom 
are coupled to each other only via interaction between their color charges, t^t^. If one will be 
able to find the limit in which this interaction becomes trivial, then the original two-dimensional 
Schrodinger equation (|3.6|) can be replaced by the system of holomorphic and antiholomorphic 
one-dimensional Schrodinger equations (|3.8| ). Moreover, this correspondence becomes exact 
in the special case of the compound states built from n = 2 and tt, = 3 Reggeons - the 
BFKL Pomeron and Odderon, respectively. For n > 3 we analyzed the Schrodinger equation 
(|3.6| ) in the multi-color limit, ^ cxd, in which holomorphic and antiholomorphic sectors 
become decoupled and the interaction occurs only between nearest Reggeons, (|3.9|) .p^ Thus, 
the expression for the energy of the compound ra— Reggeon state (|3.11|) is exact for n = 2, 3, 
but for n > 3 it contains 0{1/N'^) corrections which can be fixed from the condition of the 
5*- matrix unitarity in the 1/A^— expansion ||30[ . 



We found that holomorphic and antiholomorphic Reggeon hamiltonians (|3.1|) , ( p.2|) and 



( |3.3| ) coincide with the hamiltonian ( [4.4| ) and ( |4.5| ) of generalized XXX Heisenberg magnet for 
spin s = corresponding to the principal series representation of the S'L(2,C) group. As a 
result, the system of Schrodinger equations ( |3.8| ) turns out to be completely integrable. The 
holomorphic wave functions and the corresponding energies of the compound ra— Reggeon states 
become identical to the eigenstates and eigenvalues of the one- dimensional XXX Heisenberg 
magnet defined on the one-dimensional lattice with periodic boundary conditions and with the 
number of sites equal to the number of Reggeons. In each site of the lattice the spin operators 
are realized as generators of the conformal S'L(2,C) group acting on the Reggeon holomorphic 
and antiholomorphic coordinates. A generalized Bethe ansatz, based on the method of Baxter 
Q— operator, was developed for the diagonalization of the Reggeon hamiltonian. The spectrum 
of the n— Reggeon states, ( ^.13|) and (|5.13| ), was obtained in terms of the solution of the Baxter 

^''Another possibility would be to consider the SU (2) gauge group and take into account jlj] that for the color 
singlet n— Reggeon states — ~2/(n — l)|x) with {t'')bc — —Kabc- This means that for the SU{2) gauge 

group the Schrodinger equation ( |3.6| ) is equivalent to the system (3^) but with (anti-)holomorphic hamiltonian 
describing the long-range interaction on the one-dimensional lattice, X]fc>m H{zk, Zm)- 
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equation ( plj|) . We studied the general properties of this equation and found the relation 
between its solutions and some systems of orthogonal polynomials. For n = 2 the solutions of 
the Baxter equation were identified as series of Hahn orthogonal polynomials. These solutions 
were used to derive the spectrum of the n = 2 Reggeon state, the BFKL Pomeron. The next 
step is the solution of the Baxter for ra > 3. We might hope that the analytical solution can be 
found using the relation of Baxter equation to the orthogonal polynomials. 

Another way to solve the Baxter equation ( p.9| ) can be to explore the conformal invariance 
of the Reggeon interaction and find interpretation of the differential equation ( |6.17|) within 
the framework of two-dimensional conformal field theories. To this end, we notice the remark- 
able similarity between solution of the Baxter equation for n = 2 and the spectrum of the 
SL{2,Wi)/U{l) coset conformal field theory 



Acknowledgements 



This paper arose essentially from a joint work |T9[ with L.D. Faddeev to whom I am deeply 
indebted. I greatly benefit from stimulating discussions with J. de Boer, D. Jatkar, V.P. 
Spiridonov and G. Sterman. I am grateful for critical remarks to J. Bartels, A.V. Efremov, 
J. Ellis, D. Gross, V.E. Korepin, L.N. Lipatov, B. McCoy, A.V. Radyushkin and G. Veneziano. 
This work was supported in part by the National Science Foundation under grant PHY 9309888. 



References 

[1] S.C. Frautschi, Regge poles and S-matrix theory, New York, W.A. Benjamin, 1963; 
V. de Alfaro and T. Regge, Potential scattering, Amsterdam, North-Holland, 1965; 
P.D.B. Collins, An introduction to Regge theory and high energy physics, Cambridge Uni- 
versity Press, 1977. 

[2] HI Collaboration (T. Ahmed et al.), DESY-94-198, Nov 1994, Phys. Lett. B338 (1994) 
338; 

ZEUS Collaboration (M. Derrick, et al.) Z. Phys. C63 (1994) 391. 

[3] HI Collaboration. (K. Muller, et al.), DESY 94-112-D, Jul 1994; 

ZEUS Collaboration (M. Derrick, et al.), DESY-94-192, Oct 1994, DESY-94-143, Aug 
1994; Phys. Lett. B316 (1993) 412. 

[4] F. Low, Phys. Rev. D12 (1975) 163; 

S. Nussinov, Phys. Rev. Lett. 34 (1975) 1286; 

J. Gunion and D. Soper, Phys. Rev. D15 (1977) 2617. 

[5] A. Donnachie and P.V. Landshoff, Nucl. Phys. B244 (1984) 322; Phys. Lett. B296 (1992) 
227. 

[6] L.V. Gribov, E.M. Levin and M.G. Ryskin, Phys. Rep. 100 (1983) 1; 
E.M. Levin and M.G. Ryskin, Phys. Rep. 189 (1990) 267; 

E.M. Levin, Low x(B) physics with open eyes, FERMILAB-CONF-94-068-T, Mar 1994. 



36 



[7] G. Ingelman and P. Schlein, Phys. Lett. 152B (1985) 256; 

J. Collins et al., Measuring parton densities in the pomeron, preprint CTEQ/PUB/02, 
psu/th/136. 

L.N. Lipatov, Pomeron in quantum chromodynamics , in "Perturbative QCD", pp. 411-489, 
ed. A.H. Mueller, World Scientific, Singapore, 1989. 

[9] A.H. Mueller, Nucl. Phys. B250 (1985) 327; in QCD 20 years later, Aachen, 1992, ed. 
P.M.Zerwas and H.A.Kastrup, World Scientific, Singapore, 1993, v.l, p. 162. 

[10] M. Beneke and V.L Zakharov, Phys. Lett. B312 (1993) 340; 
V.L Zakharov, Nucl. Phys. B385 (1992) 452. 

[11] G.P. Korchemsky and G. Sterman, Nonperturbative corrections in resummed high energy 
cross sections, Stony Brook preprint ITP-SB-94-50, Nov 1994, |hep-ph/9411211| . 



[12] E. Levin, Renormalons at low x, Tel-Aviv Univ. preprint TAUP 2221-94, Dec 1994; |hep- 
I ph/9412"345| . 

[13] Handbook of Perturbative QCD, CTEQ Collaboration, ed. G. Sterman, FERMILAB-PUB- 
93/094, Apr 1993. 

[14] J. Bartels, Nucl. Phys. B175 (1980) 365. 

[15] H. Cheng and T.T. Wu, Expanding Protons: Scattering at High Energies, MIT Press, 
Cambridge, Massachusetts, 1987. 

[16] H. Cheng, J. Dickinson, C.Y. Lo and K. Olaussen, Phys. Rev. D23 (1981) 534. 

[17] E.A. Kuraev, L.N. Lipatov and V.S. Fadin, Phys. Lett. B60 (1975) 50; Sov. Phys. JETP 
44 (1976) 443; 45 (1977) 199; 

Ya.Ya. Balitskn and L.N. Lipatov, Sov. J. Nucl. Phys. 28 (1978) 822. 

[18] L.N. Lipatov, High-energy asymptotics of multicolor QCD and exactly solvable lattice mod- 
els, Padova preprint DFPD-93-TH-70, Oct 1993; JETP Lett. 59 (1994) 596. 

[19] L.D. Faddeev and G.P. Korchemsky, High energy QCD as a completely integrable model. 
Stony Brook preprint ITP-SB-94-14, Apr 1994; Phys. Lett. B342 (1994) 311. 

[20] Feynman's office: the last blackboards. Physics Today, February 1989, p. 88. 

[21] H.T. Nieh and Y.-P. Yao, Phys.Rev.Lett. 32 (1974) 1074; Phys.Rev.D13 (1976) 1082; 

B. McCoy and T.T. Wu, Phys. Rev. D12 (1975) 3257; 

C. Y. Lo and H. Cheng, Phys. Rev. D13 (1976) 1131; 
L. Tyburski, Phys.Rev. D13 (1976) 1107. 

[22] H. Verlinde and E. Verlinde, QCD at high energies and two-dimensional field theory, 
Princeton Univ. preprint, PUPT-1319, September 93. 

[23] E. Laenen, E. Levin and A.G. Shuvaev, Nucl. Phys. B419 (1994) 39. 

[24] E.K. Sklyanin, L.A. Takhtajan and L.D.Faddeev, Theor. Math. Phys. 40 (1980) 688. 



37 



[25] L.A. Takhtajan and L.D. Faddeev, Russ. Math. Survey 34 (1979) 11. 



[26] L.D. Faddeev, Stony Brook preprint, ITP-SB-94-11, Mar 1994; |hep-th/94040T^ ; The Bethe 



ansatz, Andrejewski lectures, Freie Univ. preprint, SFB-288-70, Jun 1993; in Nankai Lec- 
tures on Mathematical Physics, Integrable Systems, ed. by X.-C.Song, pp. 23-70, Singapore: 
World Scientific, 1990. 

[27] V.E. Korepin, N.M.Bogoliubov and A.G. Izergin, Quantum inverse scattering method and 
correlation functions, Cambridge Univ. Press, 1993. 

[28] J. Kwiecinski and M. Praszalowicz, Phys. Lett. B94 (1980) 413. 

[29] L.N. Lipatov, Phys. Lett. B251 (1990) 284; B309 (1993) 394. 

[30] G. Veneziano, Nucl. Phys. B74 (1974) 365; Phys. Lett. 52B (1974) 220; 
A. Schwimmer and G. Veneziano, Nucl. Phys. B81 (1974) 445; 
M. Ciafaloni, G. Marchesini and G. Veneziano, Nucl. Phys. B98 (1975) 472; 493. 

[31] V.O. Tarasov, L.A. Takhtajan and L.D. Faddeev, Theor. Math. Phys. 57 (1983) 163. 

[32] P.P. Kuhsh, N.Yu. Reshetikhin and E.K. Sklyanin, Lett. Math. Phys. 5 ( 1981) 393. 

[33] R.J. Baxter, Exactly Solved Models in Statistical Mechanics , Academic Press, London, 
1982. 

[34] E.K. Sklyanin, The quantum Toda chain, Lecture Notes in Physics (Springer) 226 (1985) 
196; Quantum Inverse Scattering Method. Selected Topics, in "Quantum Group and Quan- 
tum Integrable Systems" (Nankai Lectures in Mathematical Physics), ed. Mo-Lin Ge, Sin- 
gapore: World Scientific, 1992, pp.63-97; |hep-th/921im| . 

[35] D.P. Zhelobenko and A.I. Shtern, Representations of Lie groups (in Russian), Nauka, 
Moscow, 1983, pp.21 1-220. 

[36] E.T. Whittaker and G.N. Watson, A course of modern analysis, Cambridge Univ. Press, 
1963. 

[37] Higher transcendental functions, vols.l and 2, Bateman manuscript project, ed. A. Erdelyi, 
McGraw-Hill, 1953. 

[38] R. Askey and J. Wilson, Some basic hypergeometric orthogonal polynomials that generalize 
Jacobi polynomials , Mem. Am. Math. Soc. 319 (1985) 1-55. 

[39] V. Spiridonov and A. Zhedanov, Discrete reflectionless potentials, quantum algebras and 
q-orthogonal polynomials , preprint CRM-1928 (1993); Annals of Physics, vol. 237 (1995) 
126; Discrete Darboux transformations, discrete time Toda lattice and the Askey-Wilson 
polynomials, preprint CRM-1929 (1993). 

[40] G.P. Korchemsky, V.P. Spiridonov and A.S. Zhedanov, in preparation. 

[41] L.N. Lipatov, Talk at the 2nd Workshop on Small-x and Diffractive Physics at the Teva- 
tron, Fermilab, Sept. 1994. 



38 



[42] R. Dijkgraaf, E. Verlinde and H. Verlinde, Nucl. Phys. B371 (1992) 269; 
D. Jatkar, Nucl. Phys. B395 (1993) 167. 



39 



Figures: 





Figure 1: Unitary Feynman diagrams contributing to the structure function of deep inelastic 
scattering in the generahzed leading logarithmic approximation in the multicolor limit, iV — oo. 
Solid lines represent n reggeized gluons propagating in the t— channel. Dotted lines denote 
interaction of reggeized gluons with their nearest neighbours. Upper and lower blobs describe 
the coupling of gluons to virtual photon 7* and proton state p, respectively. For finite one 
has to add the diagrams with pair-wise interactions between n reggeized gluons. 
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Figure 2: The Bethe-Salpeter equation for the transition operator T„(c(j), describing n ^ n 
elastic scattering of reggeized gluons in the t— channel. Iterations of this equation reproduce 
the ladder diagrams of fig. |[ 
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Figure 3: The energy of the holomorphic Reggeon hamihonian for n = 3 corresponding to the 
polynomial solutions of the Baxter equation. The solid line represents the energy for n = 2 
defined in (|6.31| ). Points on this hne correspond to the degenerate solutions, gs = 0, of the 
Baxter equation for n = 3. 




Figure 4: The holomophic energy for n = 4 corresponding to the polynomial solutions of the 
Baxter equation. The solid line represents e2{h). 
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Figure 5: The holomorphic energy for n = 3 as a function of quantized q^- The maximum 
value of the energy = —6 corresponds to the degenerate solution = 0. 
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Figure 6: The holomorphic energy for n = 4 as a function of quantized ^3. The maximum 
value of energy = —8 corresponds to ga = but in contrast with fig. |^ this solution is not 
degenerate because it corresponds to q^ = 2. 
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Figure 7: The holomorphic energy for n = 4 as a function of quantized §4. The maximum 
value of the energy £4 = —8 corresponds to 54 = 2. 
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Figure 8: Quantized values of ^3 as a function of conformal weight h for n = 3 Baxter equation. 
The symmetry —>■ —q^ is in agreement with ( |6.8| ) and ( |6.9| ). The points qs = correspond 
to the degenerate solutions. The distribution of gs for n = A Baxter equation follows the same 
pattern. 
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Figure 9: Quantized values oi as a function of conformal weight hior n — A Baxter equation. 
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Figure 10: Distribution of the quantized values of versus q^ for n = 4 and h — 10. 
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